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1 Introduction 

The physical world is traditionally organized into various systems: electro- 
magnetism, perfect fluids, Klein-Gordon fields, elastic media, gravitation, 
etc. Our descriptions of these individual systems have certain features in 
common: Use of fields on a fixed space-time manifold M, a geometrical in- 
terpretation of the fields in terms of M, partial differential equations on 
these fields, an initial-value formulation for these equations. Yet beyond 
these common features there are numerous differences of detail: Some sys- 
tems of equations are linear, and some are not; some have constraints, and 
some do not; some arise from Lagrangians, and some do not; some are first- 
order, and some higher-order. Systems also differ in other respects, e.g., as 
to what fields they need as background, what interactions they permit (or 
require). It almost seems as though, in the end, every physical system has 
its own special character. 

It might be useful to have a systematic treatment of the fields and equa- 
tions that arise in the description of physical systems. Thus, there would 
be a general definition of a "field" , and a general form for a system of par- 
tial differential equations on such fields. The treatment would consist of a 
framework sufficiently broad to encompass the systems found in nature, but 
no broader. One would, for example, treat the initial-value formulation once 
and for all within this broad framework, with the formulations for individ- 
ual physical systems emerging as special cases. In a similar way, one would 
treat — within a quite general context — constraints, the geometrical charac- 
ter of physical fields, how some systems require other fields as a background. 
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how interactions operate, etc. The goal of such a treatment would be to get 
a better grip on the structural features of the partial differential equations 
of physics. Here are two examples of issues on which such a treatment might 
shed light. What, if any, is the physical basis on which the various fields 
on the manifold M are grouped into separate physical systems? Thus, for 
instance, the fields (n. Fab- P- u°') arc grouped into (n, p, u°') (a perfect fluid), 
and {Fab) (electromagnetic field). By ''physical basis", we mean in terms of 
the dynamical equations on these fields. A second issue is this: How does it 
come about that the fields of general relativity are singled out as those for 
which diffeomorphisms on M are gauge? On its face, this singling out seems 
surprising, for the diffeomorphisms act equally well on all the physical fields 
on M. 

We shall here discuss, in a general, systematic way, the structure of the 
partial differential equations describing physical systems. We take it as given 
that there is a fixed, four-dimensional manifold M of "space-time events", 
on which all the action takes place. Thus, for instance, we are not con- 
sidering discrete models. Further, physical systems are to be described by 
certain "fields" on M. These may be more general than mere tensor fields: 
Our framework will admit spinors, derivative operators, and perhaps other 
field-types not yet thought of. But we shall insist — largely for mathematical 
convenience — that the set of field-values at each space-time point be finite- 
dimensional. We shall further assume that these physical fields are subject 
to systems of partial differential equations. That is, we assume, among other 
things, that "physics is local in M" . Finally, we shall demand that these par- 
tial differential equations be first-order (i.e., involving no higher than first 
space-time derivatives of the fields), and quasilinear (i.e., linear in those first 
derivatives). That the equations be first-order is no real assumption: Higher- 
order equations can — and will — be cast into first-order form by introducing 
new auxiliary fields. (Thus, to treat the Klein-Gordon equation on scalar 
field ip, we introduce an auxiliary vector field playing the role of Vip.) It is 
my sense that this is more than a mere mathematical device: The auxiliary 
fields tend to have direct physical significance. It is not so clear what we are 
actually assuming when we demand quasilinearity. It is certainly possible to 
write down a first-order system of partial differential equations that is not 
even close to being quasilinear (e.g., (dip/dx)^ + (dip/dt)^ = ip"^). But all 
known physical systems seem to be described by quasilinear equations, and 
it is anyway hard to proceed without this demand. In any case, "first-order. 
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quasilinear" allows us to cast all the partial differential equations into a con- 
venient common form, and it is on this common form that the program is 
based. A case could be made that, at least on a fundamental level, all the 
"partial differential equations of physics" are hyperbolic — that, e.g., elliptic 
and parabolic systems arise in all cases as mere approximations of hyperbolic 
systems. Thus, Poisson's equation for the electric potential is just a facet of 
a hyperbolic system. Maxwell's equations. 

In Sect. 2, we introduce our general framework for systems of first-order, 
quasilinear partial differential equations for the description of physical sys- 
tems. The physical fields become cross-sections of an appropriate fibre bun- 
dle, and it is on these cross-sections that the differential equations are writ- 
ten. So, for instance, the coefficients in these equations become certain tensor 
fields on the bundle space. This framework, while broad in its reach, is not 
particularly useful for explicit calculations. The remaining sections describe 
various structural features of these system of partial differential equations. 
A "hyperbolization" (Sect. 3) is a casting of the system of equations (or, 
commonly, a subsystem of that system) into what is called symmetric, hy- 
perbolic form. To such a form there is applicable a general theorem on 
existence and uniqueness of solutions. This is the initial- value formulation. 
The constraints (Sect. 4) represent a certain subsystem of the full system, 
the equations of which play a dual role: providing conditions that must be 
satisfied by initial data, and leading to differential identities on the equa- 
tions themselves. The constraints are integrable if these "differential identi- 
ties" really are identities; and complete if the constraint subsystem, together 
with the subsystem involved in the hyperbolization, exhausts the full sys- 
tem of equations. The geometrical character of the physical fields has to 
do with how they "transform", i.e., with lifting diffeomorphisms on M to 
diffeomorphisms on the bundle space (Sect. 5). Combining all the systems 
of physics into one master bundle B, then the full set of equations on this 
bundle will be M-diffeomorphism invariant. This diffeomorphism invariance 
requires an appropriate adjustment in the initial- value formulation for this 
combined system. Finally, we turn (Sect. 6) to the relationships between the 
various physical fields, as reflected in their differential equations. Physical 
flelds on space-time can interact on two broad levels: dynamically (through 
their derivative-terms), and kinematically (through terms algebraic in the 
fields). Roughly speaking, two fields are part of the same physical system 
if their derivative-terms cannot be separated into individual equations; and 
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one field is a background for anotlier if the former appears algebraically in 
the derivative-terms of the latter. The kinematical (algebraic) interactions 
are the more familiar couplings between physical systems. 

It is the examples that give life to this general theory. We have assembled, 
in Appendix A, a variety of standard examples of physical systems: the fields, 
the equations, the hyperbolizations, the constraints, the background fields, 
the interactions, etc. We shall refer to this material frequently as we go 
along. Thus, this is not the standard type of appendix (to be read later, if at 
all, by those interested in technical details), but rather is an integral part of 
the general theory. Indeed, it might be well to review this material first as a 
kind of introduction. Appendix B contains a statement and an outline of the 
proof of the theorem on existence and uniqueness of solutions of symmetric, 
hyperbolic systems of partial differential equations. 

All in all, this subject forms a pleasant comingling of analysis, geometry, 
and physics. 

2 Preliminaries 

Fix, once and for all, a smooth, four-dimensional manifoldQ M. The points 
of M will be interpreted as the events of space-time, and, thus, M itself will 
be interpreted as the space-time manifold. We do not, as yet, have a metric, 
or any other geometrical structure, on M. 

We next wish to introduce various types of "fields" on M. To this end, 
let 6 A M be a smooth fibre bundle^ over M. That is, h is some smooth 
manifold (called the bundle manifold) and vr is some smooth mapping (called 
the projection mapping)] and these are such that locally (in M) h can be 
written as a product in such a way that vr is the projection onto one factor^. 
An example is the tangent bundle of M: Here, h is the eight-dimensional 
manifold of all tangent vectors at all points of M, and vr is the mapping that 
extracts, from a tangent vector at a point of M, the point of M. That the 

^We take M to be connected, paracompact, and HausdorfF. 

^See, e.g., Steenrod, The Topology of Fibre Bundles (Princeton University Press, Prince- 
ton, 1954). Note that, in contrast to what is done in this reference, we introduce no Lie 
group acting on b. 

•^This means, in more detail, that, given any point x G M, there exists an open neigh- 
borhood [/ of x, a manifold F, and a diffeomorphism C, from U x F to 7r^^[t/] such that 
-K o C, is the projection of C/ x F to its first factor. 
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local-product condition holds, in this example, is seen by expressing tangent 
vectors in terms of their components with respect to a local basis in M. 
Returning to the general case, for any point x of M, the fibre over x is the 
set of points 7r~^(a;), i.e., the set of points k G 6 such that 7r(K) = x. It follows 
from the local-product condition that each fibre is a smooth submanifold of 
h, and that all the fibres arc diffcomorphic with each other. In the example 
of the tangent bundle, for instance, the fibre over point a; G M is the set 
of all tangent vectors at x. Next, let A be any smooth submanifold of M. 
A cross-section over A is a smooth mapping A b such that tt o is the 
identity mapping on A. Thus, a cross-section assigns, to each point x of A, 
a point of the fibre over x. Typically, A will be of dimension four (i.e., an 
open subset of M) , or three. 

We interpret the fibre over x as the space of allowed physical states at 
the space-time point the space of possible field- values at x. Then 

the bundle manifold b is interpreted as the space of all field-values at all 
points of M. A cross-section over submanifold A becomes a field, defined 
at the points of A. In most, but not all, examples (Appendix A) b will 
be a tensor bundle. Thus, for electromagnetism b is the ten-dimensional 
manifold of all antisymmetric, second-rank tensors at all points of M. For 
general relativity, by contrast, b is the fifty-four-dimcnsional manifold a point 
of which is comprised of a point of M, a Lorentz-signature metric at that 
point, and a torsion-free derivative operator at that point. In both of these 
examples, the projection tt merely extracts the point of M. 

It is convenient to introduce the following notation. Denote tensors in 
M by lower-case Latin indices; and tensors in b by lower-case Greek indices. 
Then, at any point /t G 6, we may introduce mixed tensors, where Latin 
indices indicate tensor character in M at 7r(K), and Greek indices tensor 
character in b at k. For example, the derivative of the projection map is 
written (VTr)^", i.e., it sends tangent vectors in b at point k to tangent vectors 
in M at 7r(K). The derivative of a cross-section, 0, over a four-dimensional 
region of M is written (V0)a°; and we have, from the defining property of a 
cross-section, 

(V0)„"(V7r)„^ = 5a'. (1) 

A vector A" at k G 6 is called vertical if it is tangent to the fibre through 
K, i.e., if it satisfies A"(V7r)a'' — 0. Elements of the space of vertical vectors 
at K, will be denoted by primed Greek superscripts. Thus, A"' means "A is 
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a tangent vector in 6, a vector which, by the way, is vertical". Elements 
of the space dual to that of the vertical vectors will be denoted by primed 
Greek subscripts. Thus, /Iq,/ means "/i is a linear mapping from vertical 
vectors in b to the reals" . More generally, these primed indices may appear 
in mixed tensors. Note that we may freely remove primes from superscripts 
(i.e., ignore the verticality of an index), and add primes to subscripts (i.e., 
restrict the mapping from all tangent vectors to just vertical ones), but not 
the reverses. As an example of this notation, we have: (VTr)^/" = 0. 

To illustrate these ideas, consider electromagnetism. Then a typical point 
of the bundle manifold bis k = {x, Fab), where a; is a point of M and Fab is an 
antisymmetric tensor at x. A tangent vector A" in 6 at k can be represented 
as an "infinitesimal change^ in both the point a; of M and the antisymmetric 
tensor Fab'- Given such a A", the combination X°'(V7i)a"' is that tangent 
vector in M at x represented by just the "change in x"-part of A" (ignoring 
the "change in Fab" -ps^it). Such a A"' is vertical provided its "change in 
x" vanishes — so, a vertical vector is represented simply as an infinitesimal 
change in the antisymmetric tensor Fab, with x fixed. In this example, we 
might introduce the field fia'ab = fJ^a'[ab] on b, which takes any such vertical 
vector. A" , and returns, as A" fia'ab, the change in the tensor Fab at x. 

A connection on fibre bundle 6 A M is a smooth field 7a" on b satisfying 
7a"(V7r)a^ = 5a^- Givcu a connection 7^", those vectors at kg6 that can be 
written in the form ^'^7a" for some are called horizontal. Of course, there 
exist many possible connections, and so many such notions of "horizontal- 
ity". It follows directly from these definitions that, fixing a connection, every 
tangent vector in 6 at k can be written, uniquely, as the sum of a horizontal 
and a vertical vector, i.e., that every vector can be split into its horizontal 
and vertical parts. We may incorporate this observation into the notation 
by allowing ourselves the operations with primes that were previously pro- 
hibited: In the presence of a fixed connection, 7^", we may affix a prime 
to a Greek superscript (by taking the vertical projection); and, in a similar 
way, we may remove a prime from a Greek subscript. For example, we have 
7a° = 0. Note that in every case the removal and subsequent affixing of 
a prime leaves a tensor unchanged (but not so for affixing a prime and its 
subsequent removal.) 

Again consider, as an example, the case of electromagnetism. (Any other 

"^By "infinitesimal change in the point of...", we mean "tangent vector to a curve in...". 
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(nonscalar) tensor bundle would be similar.) Fix any smooth derivative oper- 
ator Va on the manifold M. Then this Va gives rise to a connection on 6, 
in the following manner. For k — {x, Fab) any point of b, and any tangent 
vector in M at x, let A" = ^"'■ja" be that tangent vector in 6 at /t represented 
as follows: "The infinitesimal change in x is that dictated by while the 
infinitesimal change in Fah is that resulting from parallel transport, via Va, 
of Fab from x along We thus specify the combination ^'^7a" for every 
and so the tensor 7a" itself. Note that we have A°(V7r)a," = which shows 
that the 7a" so defined is indeed a connection. So, the horizontal vectors 
at K, in this example are those for which "the infiniestimal change in Fab is 
exactly that resulting from parallel transport" . Clearly, every tangent vector 
in b can be written, uniquely, as the sum of a vertical vector and such a 
horizontal vector. While every derivative operator on M gives rise, as above, 
to a connection on b, there are many other connections on b (corresponding 
roughly to "non-linear parallel transport"). 

We shall not routinely make use of a connection in what follows, for two 
reasons. First, for some fields, such as the derivative operator of general 
relativity, we have no natural connection. Second, even when there is a 
natural connection (e.g., for electromagnetism) , that connection will itself be 
a dynamical variable. It is awkward having one dynamical field playing a 
crucial role in the kinematics of another. 

We now wish to write down a certain class of partial differential equations 
on cross-sections. To this end, let kA^a a-^id J a be smooth fields on b. Here, 
the index "A" lives in some, as yet unspecified, vector space. Normally, this 
vector space will be some tensor product involving tensors in M and in b, 
i.e., "A" will merely stand for some combination of Latin and Greek indices. 
But, at least in principle, this could be some newly constructed vector space 
attached to each point of b, in which case we would have to introduce a 
new fibre bundle, with base space b, to house it. Consider now the partial 
differential equation 

A;a'".(V0)^"+j^ = O, (2) 

where U — > 6 is a smooth cross-section over some open subset U of M. This 
equation is to hold at every point xEU, where k and j are evaluated "on the 
cross-section", i.e., at (f){x). Note that this is a first-order equation on the 
cross-section, linear in its first derivative. The "number of unknowns" at each 
point is the dimension of the fibre; the "number of equations" the dimension 
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of whatever is the vector space in which the index "A" hves. The coefficients 
in this equation, k^^a ci-nd Ja, are functions on the bundle manifold b, i.e., 
these coefficients may "depend on both the point of M and the field-value 
0". 

Apparently, every system of partial differential equations describing a 
physical system in space-time can be cast into the form of Eqn. (0). Var- 
ious examples are given in Appendix A. Many, such as those for a perfect 
fluid, the electromagnetic field, or the charged Dirac particle, are already 
packaged in the appropriate form. Others must be brought into this form 
by introducing auxiliary fields. In the Klein-Gordon case, for example, we 
must augment the scalar field ip by its space-time derivative, ipa, resulting 
in a bundle space with five-dimensional fibres. We then obtain, on {ip,ipa), 
a ffist-order system of equations of the form (^. For general relativity, the 
fibre over xgM consists of pairs {gab,^a), where gab is a Lorentz-signature 
metric and Va a torsion-free derivative operator at x. The curvature tensor 
arises in (|^) as the derivative of the derivative operator. Let us agree that all 
ffist-order equations on the fields (even those that follow from differentiating 
other equations) are to be included in (H). Thus, for example, Eqn. (^) is 
included for the Klein-Gordon system. Note that the only structure we are 
imposing on the physical fields at this stage is a differentiable structure, as 
carried by the manifold b. If you wish to utilize any additional features on 
these fields — e.g., the ability to add fields, to multiply them by numbers, to 
multiply them by each other, etc. — then this must be introduced, separately 
and explicitly, as additional structure on the bundle space b. For example, 
for electromagnetism, but not for a perfect fluid, each fibre has the additional 
structure of a vector space 

The present formulation of partial differential equations carries with it 
a certain gauge freedom. Let XA^b be any smooth field on b. Then Eqn. 
(ID remains invariant under adding to /cA™a the expression XA^b(^'^)a'^ and 
at the same time to Ja the expression —XA^m- That is, the solutions (p of 
(H) before these changes in k and j are precisely the same as the solutions 
after. Note that fc^'^a', (i.e., what results from contracting ^A^a only with 
vertical f ") is gauge-invariant. Furthermore, this tensor exhausts the gauge- 
invariant information, in the following sense: Given any field /ca^q satisfying 
kA^a' = kA^a': then there exists one and only one gauge transformation that 
sends kA^a to ^A^a- This gauge freedom reflects the idea that "the horizontal 
part" of the a-contraction in (^ does not really involve the derivative of the 
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cross-section, by virtue of the identity ([T|). Thus, the components of kA^a 
that participate in this part of the a-contraction are not contributing to the 
dynamics. It would be most convenient if we could somehow circumvent this 
gauge freedom, e.g., by rewriting Eqn. (0) to involve only the gauge- invariant 
part, kA^a'y of k/^a- Unfortunately, this cannot be done in any natural way 
in general. But it can be done in the presence of some fixed connection, 7^", 
on the bundle h. In fact, given a connection, we may always achieve through 
gauge a kA^^ that is "vertical in a" in the sense that it annihilates every 
horizontal vector /i". Furthermore, this requirement on kA^a completely 
exhausts the gauge freedom. Indeed, given /cA™a and connection 7^", then 
the gauge transformation with Aa™'^ = —'jb'^kA^^a: uniquely, does the job. It 
will sometimes be convenient, when a connection is available, to exploit this 
gauge-choice. 

3 Hyperbolizations 

A key feature of the partial differential equations of physics is their initial- 
value formulation, i.e., their formulation in terms of initial data and "time"- 
evolution. It turns out that this formulation can be carried out in a rather 
general setting. This is the subject of the present, and much of the following, 
section. 

Fix a partial differential equation of the form (|^), so we have in particular 
fixed smooth fields and Ja on b. By a hyperbolization of Eqn. (^, we 

mean a smooth field h'^a' on b such that i) the field h^a'kA^p' on b is sym- 
metric in a', (3'\ and ii) for each point k E b, there exists a covector 11^ in M 
at 7r(/t) such that the tensor rimh'^a'kA^is' at k is positive-definite. Note that 
the definition involves only kA^^a', and neither Ja nor the rest of k. Thus, in 
particular, the definition is gauge- invariant. Note also that the hyperboliza- 
tions at a point k G & form an open subset of a vector space. For h^^' any 
hyperbolization, and any nonzero (vertical) vector at a point, the com- 
bination f " h'^a' at that point must be nonzero. (This follows, contracting 
the positive-definite tensor in ii) with f " f ^ .) But this implies, in turn, that 
the dimension of the space of equations in (^) (that of the index "A" ) must 
be greater than or equal to the dimension of the space of unknowns (that 
of the index "a"'). So, if this dimensionality criterion fails, then there can 
be no hyperbolization. But suppose this criterion is satisfied: Can we then 
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guarantee a hyperbolization? The answer is no. In fact, there is no known, 
practical procedure, given a general partial differential equation (H), for find- 
ing its hyperbolizations, or, indeed, for even determining whether or not one 
exists. (This is essentially a little algebra problem: Given a tensor k^^a' 
at a point, what are the tensors h'^a' at that point with h^a'^A^p' symmet- 
ric?) In practice, hyperbolizations are found, in sufficiently low dimensions, 
by solving explicitly the algebraic equations inherent in i) and ii) ; and, in 
higher dimensions, by guessing. Physical considerations frequently suggest 
candidates. 

Consider again the example of electromagnetism (Appendix A). We have 
already remarked that, at point k = {x, Fab) of the bundle space b, a typ- 
ical vertical vector, which we now write Scp"^ , is represented by an infinies- 
imal change, SFab, in the electromagnetic field at x. Since the left sides of 
Maxwell's equations, (^) and (|5DD, consist of a vector and a third-rank, 
antisymmetric tensor, the index "A" lies in the eight- dimensional vector 
space of such objects. That is, a typical vector in this space can be written 

^ ^sa^^abc-^^ ^j^j^ ^abc ^ g[abc] ^ ^-^^^^ ^J^^^^ gj^^g = g > 6 = dim 

"a"', our dimensionality criterion above is satisfied.) The fields kA^a and 
are to be read off by comparing Maxwell's equations, ( pQ]) and (|30|), with the 
general partial differential equation (Q). We thus obtain 

kA"'(3'a^nrn5^^' = s'^in'SFab) + s'=^''{n[JFb,]). (3) 

Here, we have represented kA^p' by giving the scalar that results from con- 
tracting away its indices, on vectors a^, Um, and 50^'. The field Ja of (0), 
on the other hand, depends on gauge. If we choose for the gauge that de- 
termined by the derivative operator Va on M used in Maxwell's equations 
(p9D, (0), then we have Ja = 0. Now fix any vector at x, and consider the 
tensor h^a' given, in Eqn. (^TJ), as the A- index vector that results from the 
contraction . Substituting this vector for in (|), we obtain 

h^a'Sr'kA'^p'nm^^' = (5F'^™r(n^5F,,)-^t['^5F"^l(n[,5Fb,]) 

= 2en'[6F(^rSFb)m-l/^gabSFmnSF^'']. (4) 

It now follows, provided only that the vector is chosen timelike, that the 
h'^a' of (|3TD is a hyperbolization. Indeed, condition i) follows from the fact 
that the last expression in (^ is symmetric under interchange of SFab and 
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5Fab'i and condition ii) follows from the fact that, whenever rim is time- 
like with f^Um < 0, the last expression in defines a positive-definite 
quadratic form in 6 Fab- Thus, every timelike vector field t° on M gives rise 
to a hyperbolization of Maxwell's equations. In fact, this family exhausts the 
hyperbolizations in the Maxwell case. 

The situation is similar for many other physical examples. (See Appendix 
A.) Thus, the hyperbolizations of the Klein-Gordon equation are character- 
ized by two vector fields on M; and those for the perfect-fluid equation by two 
scalar fields. Even dissipative fluids^ can be described by equations admit- 
ting a hyperbolization. There are only two physical examples, as far as I am 
aware, for which there exist no hyperbolization. One is Einstein's equation, 
for which the lack of a hyperbolization is related to the diffeomorphism- 
invariance of the theory; and the other is dust. We shall return to each of 
these examples later. 

Fix a hyperbolization, h^a', of Eqn. (Q). For each point k G &, denote 
by Sk the collection of all covectors rim in M at vr(/t) such that the tensor 
nrnh'^a'^A^P' is positive-definite. Then is a nonempty (by condition ii)), 
open, convex cone. The physical interpretation of these cones will turn out 
to be the following. The tangent vectors p"" in M at tt{k,) such that p"'na > 
for all Ha G represent the "signal-propagation directions" of the physical 
field. Note that these form a closed, nonempty, convex cone at each point, 
the "dual cone" of s^- These cones depend not only on the point x of M, 
but also in general on the value of the field at x, i.e., on where we are in 
the fibre over x. In cases in which there is more than one hyperbolization, 
these cones could also depend on which hyperbolization has been selected. 
But it turns out that, for most physical examples, these cones are essentially 
independent of hyperbolization. Thus, in the case of electromagnetism, the 
signal propagation directions consist of all timelike and null vectors lying 
in one of the two halves of the light cone. In the case of a perfect fluid, the p°' 
form the "sound cone". Is it possible to isolate, via a definition, the crucial 
algebraic feature of k^^a' in such physical examples that is responsible for 
hyperbolization-independent cones? 

Suppose that, included among the various physical fields on M is a space- 
time metric, Qab- In that case, we say that the system (|]) is causal if all the 

^See, e.g., IS. Liu, I. MuUcr, T. Ruggeri, Ann Phys (NY) 169, 191 (1986); R. Geroch, 
L. Lindblom, Ann Phys (NY) 207, 394 (1991); R. Geroch, J. Math. Phys. 36, 4226 (1995). 
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signal-propagation directions are timelike or null. This is equivalent to the 
condition that each includes all timelike vectors lying within one of the 
two halves of the light cone. A perfect fluid, for example, is causal provided 
its sound speed, |^ + does not exceed the speed of light. 

Fix a hyperbolization h'^^' of Eqn. (D). This hyperbolization leads, as 
we now explain, to an initial-value formulation. By initial data we mean 
a smooth, three-dimensional submanifold S of M , together with a smooth 
cross-section, S b, over 5", such that, for every point x & S, a normal 
Um to S" at X lies in the cone Stj)a(x)- In other words, we must specify the 
physical state of the system at each point of the three-dimensional manifold 
S, in such a way that, at every point of S, all signal-propagation directions are 
transverse to 5*. Note that the role of the cross-section, (pQ, in this definition is 
to determine the cone within which the normal to S must lie. Thus, a change 
of cross-section, keeping 5* fixed, could destroy the initial-data character of 
(5, 0o). (Changing the hyperbolization could, in principle, also change the 
initial-data character, but, as we remarked earlier, it generally does not.) As 
an example of these definitions we have: If we have on M a spacetime metric 
gab with respect to which (|^) is causal, then every (S", 0o), with S spacelike, 
is initial data. 

We may now summarize the fundamental existence-uniqueness theorem 
as follows. Given initial data {S, (po), there exists, in a suitable neighborhood 

U of S, one and only one cross-section, U b, such that i) cpls = (po and ii) 

/^V[^A'"a(V0)™"+jA] = O. (5) 

Condition i) ensures that the field 0, specified over the neighborhood U of 
S, agree, on 5* itself, with the given initial conditions, 0o- Condition ii) 
ensures that the field satisfy a certain partial differential equation derived 
from (specifically, by contracting it with h^/s'). In short, the theorem 
states that we can "solve" the partial differential equation (|^), uniquely, 
subject to any given initial conditions. There is given in Appendix B a more 
detailed version of this theorem (including more information regarding the 
neighborhood U), and a sketch of the proof. This version, in particular, 
supports our interpretation of the cones in terms of signal-propagation. 

Since every solution of Eqn. (0) is automatically a solution of (^, the the- 
orem above guarantees local uniqueness of the solutions of any system, (Q), 
of partial differential equations admitting a hyperbolization. Thus, for most 
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systems of interest in physics, initial data lead to a unique local solution. Fur- 
thermore, if the hyperbolization h^a' is invertible (which holds, by the way, 
if and only if dim "A" = dim "a"'), then Eqn. (|^) is equivalent to Eqn. (j^). 
In this case, e.g., for a perfect fluid, the theorem also guarantees local exis- 
tence of solutions of @). But in many physical examples — electromagnetism 
included — h^a' is not invertible — so part of Eqn. (0) is lost in the passage 
to (^. In these cases, we cannot guarantee, directly from the theorem, local 
existence of solutions of Eqn. (H). The fate of these "lost equations" is the 
subject of the following section. 

Let us now return briefly to the example of dust. (See Appendix A.) With 
the traditional choice of fields — p (mass density) and (unit, timelike four- 
velocity) — the dust equations, (0) and (^Sj), admit no hyperbolization. This 
is perhaps surprising, for this system "obviously" has an initial-value formu- 
lation. It turns out that, if we introduce the auxiliary field wj' = VqM^, then 
the corresponding system of equations on this new set of fields, {p,u'^,wj'), 
does admit a hyperbolization. It is not clear what, if any, is the physical 
meaning of this modification. Furthermore, the hyperbolization it produces 
is apparently lost on coupling the dust with gravitation via Einstein's equa- 
tion. (This behavior is a consequence of the appearance of a Riemann tensor 
in the equations on (p, -u", wj').) What is going on physically in this example? 

4 Constraints 

Fix a partial differential equation of the form (^, so we have in particular 
fixed smooth fields kA^a and Ja on b. While much of the material of this 
section finds application to the initial-value formulation, we require at this 
stage no specific hyperbolization — nor, indeed, even the existence of one. 
A constraint at point k G 6 is a tensor c"^" at k such that 

c^('^fc^")„, = 0. (6) 

Note that the definition is gauge- invariant, and that the constraints at n form 
a vector space. A number of examples, for various physical systems, is given 
in Appendix A. For instance, the equations for a perfect fluid admit only the 
zero constraint; those for Klein-Gordon, a ten-dimensional vector space of 
constraints; and those for general relativity an eighty-four dimensional vector 
space. Maxwell's equations, on the other hand, admit a two-dimensional 
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vector space of constraints: The general c^"" is given, in this case, by Eqn. 
(|32D , where x and y are arbitrary numbers. To check that this c^" does 
indeed satisfy (^), combine it with the kA^a' for Maxwell's equations given 
by (H), to obtain 

c Ka a'Ocp = xor + ye or be- [7) 

Now symmetrize both sides over n, m. Each constraint, as we shall see, plays 
two distinct roles: It signals a differential condition that must be imposed on 
initial data for Eqn. (^, as well as a differential identity involving Eqn. (j^). 
In the case of Maxwell's equations, for example, the first role is reflected in the 
familiar "spatial constraint equations" , V • E = 0, V ■ B = 0. The second 
role is reflected in the fact that identities result from taking the divergence 
and curl, respectively, of Maxwell's equations, (^) and (|30|) . 

We begin with the first role. Fix constraint c^^. Let U b he any 
solution of Eqn. (||), defined in open U C M, and let S be any three- 
dimensional submanifold of U. Consider now the equation 

naC^1A;A™„(V0)„"+jA] =0, (8) 

at points x of S, where Ua is a normal to S* at x and the coefficients are 
evaluated at k = We first note that Eqn. (|) holds on S, for it is a 

consequence of (|^). We next claim that the left side of Eqn. (|^) involves 
only 00 = (f>\s, i-e., only restricted to 5*. To see this, first note that (po 
alone determines (V(/>)m" at points of S up to addition of a term of the form 
rimv"' . But such a term annihilates UaC^^'kA^a, by the defining equation, 
(H), for a constraint, and so does not contribute to the left side of Eqn. (||). 
What we have shown, then, is that Eqn. (^) is a "constraint equation": It 
is a differential equation on cross-sections over S that must be satisfied by 
every restriction to of a solution of Eqn. (Q). In the Maxwell case, for 
example, the two independent constraints give rise, via (P), to the vanishing 
of the divergence of the electric and magnetic fields. Note that {S, 0o) above 
need not be initial data: We have as yet introduced no hyperbolization. 

We next introduce a notion of "sufficiently many" constraints. We say 
the constraints are complete if, for any point k & b and any nonzero covector 
Un at tt{k) G M, we have 

dim(c^"n„) + dim(w"') = dim{a^). (9) 
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The first termQ is the dimension of the space of all vectors of the indicated 
form, as c"^" runs over all constraints at k. The second term is the dimension 
of the space of vertical vectors, i.e., the dimension of the fibres. The last 
term is the dimension of the space of equations represented by (|^). Eqn. 
(H) means, roughly speaking, that there are at least as many equations as 
unknowns in Eqn. (Q), and that any excess is taken up entirely by constraint 
equations, (§). This interpretation will be made more precise shortly. 

The constraints are complete for the vast majority of physical examples. 
(See Appendix A.) Thus, Eqn. (P) reads, for the perfect-fiuid equations, 
"0 + 5 = 5" ; for Maxwell's equations, "2 + 6 = 8" ; for the Klein-Gordon 
equations, "6 + 5 = 11". For Einstein's equation, the constraints are not 
complete: Eqn. (|^) reads "64 + 50 = 110". This, as we shall see later, is 
related to the diffeomorphism-invariance of the theoryQ. Is there some simple 
characterization of those tensors k^^^' that yield complete constraints? 

The second role of a constraint is in signaling a differential identity in- 
volving Eqn. (^. The idea here is very simple. Eqn. (H) is to hold at every 
point X of some open subset U of M. Taking the x- derivative, Vn, of this 
equation, and contracting with any constraint c"^", we obtain an equation in- 
volving the first and second derivatives of the cross-section. But, as it turns 
out, the second- derivative term drops out, by virtue of @, and so we are 
left with an algebraic — in fact, quadratic — equation in the first derivative, 
(V0)m°, of the cross-section. That is, we obtain an integrability condition 
for Eqn. @. If this integrability condition holds as an algebraic consequence 
of Eqn. (|^), then we say our constraint is integrable. 

Unfortunately, all this becomes somewhat more complicated when writ- 
ten out explicitly. Fix any (torsion- free) derivative operator, Vq,, on the 
manifold b, such that the derivative of every vertical vector field is verti- 
cal. (Such always exists, at least locally, by the local-product character of 
the fibre bundle.) Extend^ this operator to mixed fields on b by demanding 

^Note that this term can be — and in examples (such as Klein-Gordon) frequently is — 
less than the dimension of the vector space of constraints. 

^We remark that there exist examples (though apparently no physically interesting 
ones) of a tensor fc^i™^ admitting a hyperbolization, but whose constraints are not 
complete. 

^This is done as follows. Any field ^™ can be written in the form (Vtt)/?™^'^, 
uniquely but for the freedom to add to any vertical vector field. Now define 
Vq,^™ — iy ir) p™''^ aS,^ 1 noting that the right side is invariant under this freedom. Note 
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V/3(V7r)a'" = 0. Then the operator "derivative along the cross-section" is 
(V0)„"Vq,. Applying this operator to (|^), and contracting with any con- 
straint c"^", we obtain 



c 



^"(V^fc^"J(V0)/(V0)„" + c^"(V^ja)(V</))/ = 0. (10) 



In the derivation of Eqn. ([ToD , there arises initially the term [c^^fc^'' 



[(V0)„^V^(V0)m"], involving the second derivative of the cross-section. To 
see that this term vanishes, first note that the index "a" in the second factor 
is vertical (contracting with (VTr)^'*), and so only the antisymmetrization 
of this factor over n,m contributes (by definition of a constraint), yielding 
zero (by the torsion-free character of V^). Eqn. (|10D is our integrability 
condition. We say that constraint c^"" is integrable if Eqn. (|TOD is an algebraic 
consequence of Eqn. (H). What this means, in more detail, is that the left 
side of Eqn. (^) is some multiple of the left side of (Q) plus some multiple of 
the difference between the two sides of the identity (|l]), where each of these 
two multiplying factors is an expression linear in (V0)a"- Writing this out 
and equating powers of (V0)a", we conclude: Constraint c^"' is integrable if 
and only if there exist tensors <J^"^a and cr'^b"^a, with a"'a"^a = 0, such that 

V,(c^"fcA"/3)+V;3(c^"^A"«) 

= J^"^ + + a'"/^(V7r),^ + a",'",(V7r)/, (11) 

where we have set kA^a = kA^a + |jyi(V7r)a™. Applying a prime to both 
"a" and in this equation, and using (|), we obtain 

2V[„,(c^'"A;m"^,]) = a^™„,fc^V + (12) 



This part of ( [111) is manifestly gauge-invariant (involving only kA^a', and not 
Ja or the rest of k), and independent of the derivative operator Vq, (involving 
only the "vertical curl"). What remains of Eqn. ([iTD is essentially one scalar 
relation, expressing the divergence of c^^'Ja in terms of other fields. Is there 
some simple way to write this remaining relation, e.g., a way that separates 
its physical content from the gauge freedom inherent in {k, j), Vo, and the 
cr's? In electromagnetism, to take one example, Eqn. (|l^) is satisfied with 
a^"^a' = 0. What remains of Eqn. (|TT|) in this example is just the vanishing 
of the divergence of the electric charge-current. 

that this extension of Vq to fields with Latin indices is unique. 
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Failure of integrability would mean that we have somehow failed to in- 
clude in all the relevant conditions on the first derivative of the cross- 
section. The standard procedure, in such circumstances, is, first, to enlarge 
Eqn. (0) to encompass the additional conditions on (V0)m°. Then look for 
any additional constraints arising from this enlargement, and if any of these 
fail to be integrable, enlarge Eqn. (Q) further, etc. Unfortunately, it is not 
clear, in the present general context, how to implement this procedure. How 
do you "enlarge" a system, (Q), linear in (V^)™"? to encompass a quadratic 
relation (|10D? 

Nowhere in this section so far have we introduced a hyperbolization. It is 
perhaps striking that so much of the subject of constraints can be carried out 
at this level, for it is largely in their interaction with hyperbolizations that 
constraints come to the fore. We turn now to this interaction. Fix, therefore, 
a hyperbolization, h'^a'i for Eqn. (|^). 

Let c"^" be a constraint. Then Eqn. (Q) holds for the restriction, 00 = (/>|s, 
of any solution, 0o, of Eqn. @ to any three-dimensional submanifold, 5, of 
M. So, in particular, this equation holds when (5,00) are initial data, i.e., 
when the normal to 5* at each point k lies in the cone s^f^i^^)- Thus, given 
initial data, (S", 0o); "we have no hope of finding a corresponding solution of 
Eqn. (H) unless those data satisfy Eqn. (P) for every constraint c"^". Eqns. 

become constraint equations on initial data. 

Next, fix K G 6 and G s^. Then, we claim, for any constraint c"^" and 
any (vertical) vector , we can have UaC^"" = f " h^a' only if each side is zero. 
Indeed, this equality implies {naC^"')kA^f3'v'^ rim = ("^^ h^a')kA^ p'V^ '^m- But 
the left side vanishes (by definition of a constraint), while vanishing of the 
right side implies = (by G s^). What we have shown, in other words, 
is that the subspace of vectors of the form naC^"" with c^"" a constraint, and 
that of vectors of the form f " h^a' with f " vertical, have only the zero vector 
in common. But this implies that the left side of Eqn. (^ is less than or equal 
to the right side. That is, in the presence of a hyperbolization, "half of Eqn. 
(P) is automatic. Now suppose that the constraints are complete, i.e., that the 
full equality (^) holds. It then follows that our two independent subspaces — 
{naC"^""} and {v^ h'^a'} — in fact span the space of all vectors of the form a"^. 
What this means, in geometrical terms, is that the "constraint components" 
of Eqn. (^ — the results of contracting it with vectors of the form UaC'^"' — 
and the "dynamical components" of Eqn. @) — the result of contracting it 
with h^a' — together comprise the whole of Eqn. (^. Completeness, in the 
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presence of a hyperbolization, means the absence of any "stray equations" in 

Finally, we return to the issue, raised in the previous section, of when 
there is an initial- value formulation for the full equation (^. Fix a hyper- 
bolization h^a' for this equation, and suppose that its constraints are both 
complete and integrable. Let (5, ^o) be initial data, and suppose that these 
data satisfy all the constraint equations of the type (H) (for if not, then there 
is certainly no evolution of these data). By the general existence- uniqueness 

theorem (Sect. 3 and Appendix B), there exists a solution, U b, of the 
evolution Eqn. (||), with 01^ = (po, where U is an appropriate neighborhood 
of the three-dimensional submanifold 5* of M. We now claim that, under 
certain conditions, this cross-section (p satisfies our full system, (0), of par- 
tial differential equations. It is convenient, for purposes of this paragraph, to 
introduce upper-case Greek indices to lie in the vector space of constraints; 
so, in this notation, we have a single constraint tensor, c^^'r- Denote the left 
side of Eqn. (|^), evaluated on the cross-section 0, by I a- Thus, we have 
h^a'^A = everywhere in U (by (^), and Ia = on S (by completeness); 
and we wish to show = everywhere in U. To this end, consider the 
expression 

(c^'^r + cr'""'r/^^«') V„^/a, (13) 

where a^"'r is any field on b. We claim that this expression is, everywhere 
in U, a multiple of I a- Indeed, the first term in parentheses leads to such 
a multiple since, by integrability, C^^rVm-^A is a multiple of I a', and the 
second term also leads to such a multiple, using h^^a'^A = and differen- 
tiating by parts. To summarize, we have shown so far that Ia vanishes on 
S, and satisfies a certain first-order, quasilinear (in fact, linear) partial dif- 
ferential equation arising from the expression (pIS]). Since this differential 
equation clearly has as one solution Ia = 0, we can conclude that = in 
a neighborhood of S if we can show local uniqueness of its solutions. 

The most direct way to prove local uniquenss of solutions of a partial dif- 
ferential equation is to show that it admits a hyperbolization. In the present 
instance, tensor h^^ is a hyperbolization of the differential equation resulting 
from (|I3]) — for some choice of o"*"" a — provided h^^ has the following prop- 
erty: The expression h^^c^"^rVA'WB, for all v,w with VAh^a' = WAh^a' = 0, 
is symmetric under interchange of v and w, and, contracted with some n^, is 
positive-definite. When — in terms of the original kA^a' and its hyperboliza- 
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tion h'^a' — does such a hyperbolization h"^^ exist? In physical examples — 
e.g., in electromagnetism — such an h^^ does indeed exist, and so we have in 
these examples uniqueness of solutions of the equation resulting from (]TB|), 
and so an initial- value formulation for the full system (^. Is there any simple, 
reasonably general, condition on kA^a'^ h^a' that will guarantee existence of 
a hyperbolization h^^7 Are there interesting cases in which uniqueness of so- 
lutions of the differential equation arising from (13) must be shown by some 
other method? 



5 The Combined System — Diffeomorphisms 

In the preceeding sections, we have been analyzing the structure of the partial 
differential equation describing a single physical system. This analysis was to 
be applied separately to the electromagnetic field, a perfect fluid, or whatever. 
However, in the real world, all these systems coexist on M, normally in 
interaction with each other. We now consider the system that results from 
combining all these subsystems. 

Once again, we have a smooth fibre bundle, B M, over the four- 
dimensional space-time manifold M. Now, however, the fibre in B over 
x e M represents the possible values at x of all possible physical fields in 
the universe. Thus, this fibre would include an antisymmetric tensor (for 
electromagnetism), a Lorentz metric and derivative operator (for gravity), 
two scalar and one vector field (for a perfect fluid), etc. Note that we are 
implicitly assuming that the space that results from combining these fields 
is finite-dimensional, and that we are somehow capable of "finding" it. One 
or both of these assumptions may be incorrect. In any case, we imagine 
that we have constructed such a bundle. Again, a cross-section, M ^ B, 
of B represents an assignment of a complete physical state (of everything) 
to each point of space-time, i.e., a statement of the entire dynamics of the 
universe. And, again, we impose on such cross-sections the general first-order, 
quasilinear partial differential equation, 

i^A'"a(V<l')^" + Ja = 0. (14) 

In the analogous equation for a single system, (^, we allowed the coeffi- 
cients, k^^a and Ja, to be arbitrary (smooth) fields on the bundle manifold 
b. That is, we allowed these coefficients to depend on both "the point of the 
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space-time manifold M and the value of the field assigned to that point". 
But, in the context of this combined system, an explicit dependence of the 
coefficients, Ka^u and J^, on the point of M is, we suggest, inappropriate. 
After all, we identify the points of the manifold M, not by somehow "per- 
ceiving them directly" , but rather more indirectly, by observing the various 
physical fields on M. So, for instance, the physical distinction between two 
points, X and of M rests on the difference between the values of some 
physical field at x and at y. (This issue did not arise in the context of a 
single system, for there x-dependence of k and j could arise through other 
physical fields, not included in the dynamics of (g).) In any case, we expect 
that the coefficients of K and J in ([l^ will depend explicitly only on the 
fibres of B, with any dependence on the point of M arising only implicitly 
through the cross-section. Unfortunately, this expectation — at least, as it is 
stated above — does not make mathematical sense! The problem is that our 
fibre bundles are not naturally products, and so there is no such thing as "a 
function only of the fibre- variables, independent of the base-space variables" . 
We must therefore proceed in a different way. 

We now demand that, as part of the physical content of the bundle B, 
there be given on it the following additional structure. To each diffeomor- 
phism V on the manifold M, there is to be assigned^ a lifting of it to a 
diffeomorphism, T>, on the manifold B. By "lifting", we mean that T> must 
satisfy 11 o D = D o H, i.e., that V must take entire fibres to fibres, such 
that the induced diffeomorphism on M is precisely the original V. We fur- 
ther require of these liftings that they respect the group structure of the 
M-diffeomorphisms, i.e., that (id^) = (ids) and {T> o T)') = T>oT>'. In short, 
we must specify how the physical fields "transform" under diffeomorphisms 
on M. In the examples of Appendix A, the fibres consist of tensors, spinors, 
derivative operators, etc., and on such geometrical objects there is a natural 
action of M-diffeomorphisms. Indeed, we claim that it is this "transforma- 
tion behavior" that endows such fields with a geometrical content in terms of 
M. For instance, given a point of the tangent bundle of M, the direction in 
M in which the vector "points" can be read out from the M-diffeomorphisms 

^This demand rules out most gauge theories other than electromagnetism (Appendix 
A). What happens in these theories is that to each diffeomorphism on M is assigned 
a number of hftings to B. Indeed, the collection of all liftings assigned to the identity 
diffeomorphism on M is called the gauge group. Presumably, much of what follows could 
be generalized to include such gauge theories. 
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whose lifts leave this point invariant. Consider now the bundle B that results 
from combining all the examples of Appendix A. Lift diffeomorphisms from 
M to this B by combining these liftings for all the individual examples. 

We now have the machinery to express the idea that the coefficients in 
Eqn. (|T^ be "functions only of the physical fields". We demand that, for 
every diffeomorphism P on M, its lifting T> leave K^^a and Ja invariant 
(noting that this makes sense for fields having indices in both B and M), up 
to gauge. It then follows that, for $ any cross-section satisfying Eqn. (|1^) 
and V any diffeomorphism on M, the "transformed cross-section", Vo^oV, 
is also a solution. 

The "infinitesimal version" (with diffeomorphisms replaced by vector fields) 
of all this is the following. To each (smooth) vector field on M, there is 
to be assigned a lifting of it to a vector field, on B. By "lifting" we mean 
that ^"(Vn)Q," = C,"'- We require of these liftings that they be linear (i.e., 
that {c^ + rj) = cC, +17, for c constant), and Lie-bracket preserving (i.e., that 
[^,77] = [^,77]). Invariance of the coefficients in ([T^) under these infinitesimal 
diffeomorphisms now become^ 

^^^A^a = ^A^bC^'n)^, ^^Ja = -^A^m, (15) 

for some field A^^b on B, and for every vector field on M. We shall 
further assume that results from through the action of some differential 
operator[^ 

r = <5"''"^-'"\v^,---v^,r + ---. (16) 

In (|TBp, we have written out only the highest-order term. Its coefficient, 
^ami-ms^ _ {mi-ms) ^ ^ jg gQj-Qg smooth field ou B, independent of the 

derivative operator employed in (^). Note that the index "a" of 6 is vertical, 
as follows from the definition of a lifting. As an example, consider the system 
resulting from combining all the examples of Appendix A. In this case the 
highest order appearing in the expression, (0), for is s = 2, and this 
order occurs only for the derivative operator of general relativity. A typical 
vertical vector in the bundle of derivative operators is given by 50° = Sr°-bc 

^°The infinitesimal version of the transformation of solutions of (|l^) under diffeomor- 
phisms becomes that, for every its lifting be a linearized solution of (p^. 

^^It is possible that this assumption follows already from the general properties above 
of the liftings of diffeomorphisms. 
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(Appendix A). Then the 6 of (|T6|) becomes 

5"'"^™% = (17) 

reflecting the action of a Lie derivative on a derivative operator. 

The equation, (p!^, for the combined system can never admit any hyper- 
bohzation. To see this, first note that, given any three-dimensional submani- 
fold S of M, there always exists a diffeomorphism "D on M that is the identity 
in an arbitrarily small neighborhood of S, but not outside that neighborhood. 
But now, were there a hyperbolization, then the transformation of solutions 
$ of ( [I^ ) by such a diffeomorphism would violate the uniqueness theorem for 
hyperbolic systems (Sect 3 and Appendix B). Another way to see this is to 
note that Eqns. (|I5|) and (^) together imply 

^^(m^^la'lmi-™.)^, = Q, (18) 

which implies in turn that, for any n^, the tensor KA^a'^m is annihilated 
on contraction with 5° ™^"'™*r^mi ■ ■ '"^ms- But a nonzero vertical vector an- 
nihilating KA^a'^m precludes a hyperbolization. 

How, in light of this observation, are we ever to recover an initial-value 
formulation in physics (i.e., in (|T^) )? The answer is that we must formulate 
a modified version of the initial-value formulation, in which, given suitable 
initial data, solutions of ([T^) will always exist, but will be unique only up 
to the diffeomorphism freedom. We now describe one scheme (suggested by 
general relativity) to implement this program. There may well be others. 
The idea is to supplement Eqn. (14) on the cross-section $ with a second 



system of equations, of the form 

i^aa(V$)," + Pab = 0, (19) 

where and Vab are smooth fields on B. We wish to so arrange matters 
that i) Eqns. (|T^) and (|T^), taken together on B, have an initial- value 
formulation, and ii) Eqn. (|19D can always be achieved, in a suitable sense, 
via the diffeomorphism freedom. 

We first consider the issue of an initial- value formulation for Eqns. ([T^), 
([lyp. A hyperbolization for the system (|T^), (|T^ consists of fields H^a' and 
/q,/"^" on B such that the expression 

H^o.'Ka'^p' + /.'"'V,^, (20) 
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is symmetric in indices a', and is positive-definite on contraction with 
some covector Um at each point. What of the constraints? Let us demand 
that Vaa have rank four (i.e., that ^"'I'aoi = imphes = 0). Then (|I^) 
represents a sixteen-dimensional vector space of additional equations on 
The general constraint for Eqn. (|19D is given by c"^" = x"^", with x"*" = x"'''"]. 
Thus, the dimension of the vector space of constraints is twenty-four, while, 
for any fixed n„, the dimension of the space of all ^^"n^ is just twelve. Now 
consider Eqn. (^), the condition for completeness of the constraints. As we 
have just seen, the act of supplementing Eqn. (|lj) by (|l^) increases the first 
term in Eqn. by twelve, the second term by zero, and the third term by 
sixteen. Thus, in order that the system ([I^) , (p!9|) be complete, the original 
system, must have yielded a left side of Eqn. (|^) exceeding the right 

side by exactly four. Finally, integrability of the constraints arising from Eqn. 
(|TUp yields a system of partial differential equations on the coefficients 
and Pah- The geometrical content of these equations is that ( ]T9| ) be precisely 
the requirement that the cross-section $ lie within a certain submanifold V 
of B. This V has codimension four (i.e., dimension four less than that of B), 
and its tangent vectors are those ^" with ^"t'ao = 0. 

We may summarize the discussion above as follows. Let there be given 
a field Vaa' on B satisfying the following three conditions: i) for some fields 
H and /, the expression ( pO] ) is symmetric, and, contracted with some n^, 
is positive-definite; ii) at each point of -B, has rank four, so the space 
of vectors f " at each point with v"" Vaa' = has codimension four; and iii) 
these vector spaces can be integrated to give submanifolds, of codimension 
four, within each fibre. Given such a field Vaa') choose any submanifold V 
of B, of codimension four, such that V intersects each fibre of B in one 
of the submanifolds in iii) above. Now restrict the cross-section $ to this 
submanifold V , yielding an equation of the form (p^. By construction, the 
system (0), ( |19|) admits an initial- value formulation. 

We next turn to the issue of "achieving" Eqn. (|T^) via diffeomorphisms. 
Let $ be any solution of Eqn. (0), and V any submanifold of 5, as described 
above, such that, over some three-submanifold S of M, the cross-section 
$ lies within V . We wish to find a diffeomorphism on M that sends the 
entire cross-section $ to lie within the submanifold V . We may choose our 
diffeomorphism, together with its first (s — 1) derivatives, to be the identity 
on S", but it must begin to differ from the identity on evolution off S. We 
shall be able continually to adjust $, via a diffeomorphism, to lie within 
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V provided we can generate, via (|16|), any transverse to V. But this is 
precisely the statement that the operator 

,,5-'---%v^,---v„,r (21) 



be hyperbohc. By "hyperbohc" here, we mean that the system that results 
from introducing as auxiliary fields the first (s — 1) derivatives of admits 
a hyperbolization in the sense of Sect. 3. See the discussion at the end of 
Appendix A. 

To summarize, we may recover an initial-value formulation^ for Eqn. 
( p!^ provided we can find a field Uaa' on B satisfying conditions i)-iii) above. 



together with: iv) the operator in (|2T|) is hyperbolic. As an example (in fact, 



the example) of this scheme, consider Einstein's equation (for the derivative 
operator) in general relativity^ In this example Uaa' is given by 



l^aa'Sr = gabg^'^ST",,. (22) 

This ly indeed satisfies the four conditions listed above. First note that, for 
Einstein's equation, the left side of (H) does indeed exceed the right side by 
four. The H^'^ of Eqn. (^) is given by 



_l,abcd^Y^^agrst'), (23) 

for any timelike, and s'^'"^'^ = s^^*^)*^"^) positive-definite in its two index pairs. 
Using ( p!7D and (p2D, the operator of (^) is just the wave operator. Note that 
the field Uaa' of (^) is diffeomorphism- invariant. Thus, the breaking of the 
diffeomorphism invariance takes place solely through the choice of Uaa, i-e., 
the choice of submanifold V. Are there any other t'aa''s that work for general 
relativity? 

^^There is a somewhat more elegant, if less accessible, way to formulate this. Consider 
s = 2 in Modify the bundle B to include, in each fibre, a copy of the twenty- 

dimensional manifold of all tensors pa'' at all points of M. Think of a cross-section of 
this modified bundle as including a diffeomorphism on M and a "candidate" (pa^) for the 
derivative of this diffeomorphism. Eqns. (|lj) and ([l^) can be combined as a hyperbolic 
system on this modified bundle. The "unfolding of the diffeomorphism" is then already 
incorporated in the modified bundle. 

iSpQj. recent work on the initial-value formulation of Einstein's equation, see, e.g., S. 
Fritelh, O. Reula, Comm. Math. Phys. 266, 221 (1994); Y. Choquet-Bruhat, J. W. York, 
CR Acad Sci, Paris (to appear). 
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The diffeomorphisms, of course, act simultaneously on all the fields in the 
combined bundle B. What, then, is the feature that singles out the field Va 
of general relativity (as opposed, say, to the field Fab of electromagnetism) to 
be the one for which the diffeomorphisms are taken to be the "gauge"? The 
answer, we suggest, is the following. Let be any covector in M which, 
contracted into the expression (pO]), yields a positive-definite quadratic form. 
Now apply this quadratic form to the vertical vector 



V 



°' = 5"'"^-"^^ri„,---n„,r- (24) 



Then the first term arising from ( pOD vanishes, by Eqn. p8|), and so, by 
positive-definiteness, we must have I'aa'v"' nonzero. We may restate this 
observation as follows: That field whose diffeomorphism-behavior involves 
the highest number of derivatives of in Eqn. ([T6|) (i.e., that is repre- 
sented by vertical vectors of the form (p^) is also the field restricted by 
the gauge- fixing equation, (|19D (i.e., that corresponds to vertical vectors not 
annihilating Uaa')- In the case of the combined bundle B resulting from the 
examples of Appendix A, the highest number of derivatives of arising from 
diffeomorphism-behavior is s = 2, and the field having this behavior is, of 
course, the derivative operator Va- In this manner, the derivative operator of 
general relativity acquires the diffeomorphisms as its gauge group. These re- 
marks have the following curious consequence. Suppose that, at some time in 
the future, there were introduced a new physical field, having order s = 3 in 
Eqn. (|16|). Then, apparently, that new field would take the diffeomorphisms 
as its gauge group; leaving no "gauge freedom" in general relativity. 



6 Physical Systems — Interactions 

We adopt the view that the combined bundle B, with its partial differential 
equation ( P^^D and action of the diffeomorphism group, comprises all there 
is in the (nonquantum) physical world. But, by contrast, we do not view 
our world as such a single entity. Rather, it appears to be divided into 
various "physical systems". For example, one such system is comprised of 
the electromagnetic field Fab alone; another, of the four- velocity u"', mass 
density p, and particle-number density n for a perfect fiuid. We do not, 
e.g., organize these four fields into one system {Fab,n), and another {u°',p). 
We then think of these individual systems as "interacting" with each other. 
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Interactions take place on a number different of levels. Consider, for example, 
the electromagnetic field Fab- In the absence of a space-time metric, there is 
no natural choice for the electromagnetic field tensor {Fab, or F*^*, or some 
density?); and, even if some such choice were made, there is no way to write 
down Maxwell's equations. (Note that neither of these assertions is true with 
the roles of the metric and electromagnetic field reversed.) This is one level 
of interaction. On a different level is the interaction of a charged fluid on the 
electromagnetic field through the appearance of the the fluid charge-current 
in Maxwell's equations. These structural features of the world — the notion 
of physical systems and their various levels of interaction — must somehow be 
"derived" from Eqn (|14D on B — at least, if our view of the primacy of (|l^) 
is to be maintained. How all this comes about is the subject of this section. 

Let 6 A M be a flbre bundle. By a quotient bundle of b, we mean a 
smooth manifold b, together with smooth mappings b b ^ M, such that 

i) TT o C = 7r, and ii) 6 A M is a flbre bundle (over M), and 6 6 is a flbre 
bundle (over b). Thus, a quotient "inserts a manifold b between b and M, 
in such a way that there is created a flbre bundle on each side of 6". The 
following example will illustrate both the mathematical structure and the 
types of applications we have in mind. Let b be the bundle whose flbre over 
X G M consists of pairs {gab, Fab), where gab is a Lorentz-signature metric 
at X, and Fab an antisymmetric tensor (the electromagnetic fleld). Now let 
6 A M be the bundle whose flbres include only the Lorentz metrics, and let 

b b he the mapping that "ignores Fab" ■ This b ^ b ^ M, we claim, is a 
quotient bundle. Furthermore, it reflects the natural relationship between F 
and g, i.e., that it is meaningful to discard Fab while retaining gab, but not 
to discard gab while retaining Fab- 
Let 6 6 A M be a quotient bundle. Fix any cross-section, M b, of 
b. We now construct, using this 0, a new flbre bundle, b A M, as follows. 
For the manifold b, we take the submanifold of b (i.e., the set of 

points of b lying above the flxed cross-section of 6); and for the mapping 
TT we take the restriction to the submanifold b of the projection vr. The 
bundles b ^ M and b ^ M represent a kind of "splitting" of the original 
bundle b M. For instance, we have (dim flbre b) + (dim flbre b) = (dim 
flbre b) . Every cross-section, 0, of b yields both a cross-section of b (namely, 
= C o 0), and a cross-section of b (namely, 0). And, conversely, cross- 
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sections of b and b combine to form a cross-section of b. But — and this is the 
key point of the construction — the bundle b requires for its very existence 
a given cross-section of b. We illustrate this construction with our earlier 
example (with 6-fibres {gab, Fab), and 6-fibres {gab))- Fix a cross-section of b 
i.e., fix a metric field cjab on M. Then the submanifold b of b consists of all 
{x, (jab. Fab)- That is, the metric at each x G M is required to be the fixed 
(jab there, while the electromagnetic field Fab at x remains arbitrary. So, b in 
this example is the bundle over M of electromagnetic fields, in the presence 
of the fixed background metric gab- Clearly, a cross-section of the bundle b 
(metric field gab), together with a cross-section of b (an electromagnetic field 
in the presence of background metric (jab), yield a cross-section of the original 
bundle b; and conversely. 

The notion of a quotient bundle captures the idea of one set of physical 
fields serving as the kinematical background for another. Thus, the metric 
serves as the kinematical background field for the electromagnetic field; and 
the metric and electromagnetic fields together serve as the kinematical back- 
ground for a charged scalar field. We turn now from kinematics to dynamics. 

Fix a fibre bundle b M, and a quotient bundle thereof, b b M. 
Then, as we have just seen, the cross-sections of b "split", in the sense that 
specifying a cross-section of 6 is equivalent to specifying a cross-section 
of b, together with a cross-section of the bundle b derived from 0. Next, 
let there be specified a system of quasilinear, first-order partial differential 
equations, (0), on cross-sections of b. When does this equation also split 
into separate equations on the cross-sections, and 0, that comprise 0? We 
are here concerned only with the first term in (0) — the dynamical part of the 
differential equation. The remainder will be discussed shortly. So, fix a field 
kA^a' on b: We wish to split it into corresponding fields on b and b. There is 
a natural choice for the field on 6, namely the restriction to the submanifold 
6 of 6 of the given field k^^a' on b. On taking this restriction, only some of 
the components of k^^a' survive, namely, those represented by contraction 
of kA^a' with vectors f ° tangent to the submanifold b. The remaining 
components — those lost under this restriction — must now be recovered from 
a suitable field, k£^s,i on b. This recovery will occur provided i) the pullback 
of kj^"^a' from 6 to 6 is a linear combination of k^^a' on b, and ii) this 
pullback, together with the restriction of kA^a' to b, exhausts kA^a'- The 
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first condition means, in more detail, that, for some field /i^"^ on b, we have 

(VC)a'"'A;i"<i'=^/fcA'"a', (25) 

where the right side is evaluated at k G &, the left side at C(^) ^ b. The left 
side of (|25|) is the pullback of k^"^^> from b to b via the mapping (; the right 
side, some linear combination (with coefficients fi) of kjs^a'- This condition 
guarantees that the dynamical part of the differential equation to be imposed 
on b will come from the dynamical part of the differential equation originally 
given on b. The second condition means, in more detail, that, for any vector 
such that the restriction of a'^kA^a' to b vanishes, we have = r^fj,^"^ 
for some t^. In other words, what is lost on restriction (to b) must be 
regained via the pullback (from b). Note that such a A; on 6, if it exists, is 
unique. 

Given fibre bundle b A M, and field kA^a' on fe, by a reduction of this we 
mean a quotient bundle 6, with field k^a' on 6, satisfying the two conditions 
above (in the case of ii), for every cross-section 0). To illustrate this defini- 
tion, we return again to our earlier example. On the metric-electromagnetic 
bundle b above, introduce the field kk^a' arising from the metric-Maxwell 
equations: 

Va^?fec = 0, V'F,b = 0, V[aFfe,] = 0. (26) 

Let i) be the quotient bundle above (in which only the metric gab is retained). 
On 6, introduce the field kj^a' arising from the first equation in (^Bj). (Here, 
we are making essential use of the fact that no electromagnetic field appears 
in this equation.) This 6, k^a'-, "we claim, is a reduction. Indeed, condition 
i) follows from the fact that, whenever cross-section {gab, Fab) of the bundle 
b satisfies the full system (^), then the corresponding cross-section gab of b 
satisfies the first of these equations. Condition ii) follows from the fact that, 
given a field gab (cross-section of b) satisfying the first equation, and then an 
Fab (cross-section of b) satisfying, in the presence of that gab as background, 
the last two equations, then the full set, {gab, Fab), satisfies the full system 
(p6[). Note, by contrast, that there is no reduction with the roles of the 
electromagnetic field and the metric reversed. 

Let us now return to the combined bundle B, with its partial differential 
equation (^). This {B, Ka^o') will, of course, have numerous reductions; 
and the {B, Kj^a') that result may have further reductions; and so on. Any 
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fields lost (i.e., incorporated in bundle B) in such a reduction will be called 
a physical system; while the fields remaining (in the bundle B) will be called 
background fields for that physical system. These definitions, we suggest, cap- 
ture the way in which fields are grouped together in physics, and the sense in 
which the equations for some fields require as prerequisites other fields. For 
example, the fields (m", p, n) for a perfect fiuid form a physical system (but 
not, e.g., (p, n) alone); with background the space-time metric. The charged 
Dirac field is a physical system, with background the electromagnetic and 
metric fields together; and the electromagnetic field is a physical system, 
with background the metric field. This definition also produces a couple of 
minor surprises. For the Klein-Gordon equation, the scalar and vector fields, 
ip and ipa, form separate physical systems, the former having no background 
fields; the latter, just the metric field. (Thus, neither of these has the other 
as background!) Similarly, the metric and derivative operator of general rel- 
ativity form separate physical systems, the metric having no background, 
the derivative operator the metric as background. Note, from the examples 
of Appendix A, that the metric is a background for virtually every physical 
system (the sole exceptions being those, such as the Klein-Gordon -0, having 
such a wide variety of hyperbolizations that every covector in M lies in the 
cone Sk for at least one of them). This feature presumably refiects the fact 
that, in order that there be a hyperbolization for the combined system (0), 
it is not enough merely that there be a hyperbolization for each individual 
physical system making up B. These individual hyperbolizations must also 
be such that they have in their various a common n^. In order to achieve 
this, the individual physical systems must somehow arrange to "communi- 
cate" with each other what their dynamics is. That communication takes 
place by sharing the space-time metric gab as a background field. 

So far, we have focussed exclusively on the "dynamical part" of Eqn. 
( p!4D — the K^^a'- We turn now to the remainder of this equation — the Ja- 
Roughly speaking, whenever that portion of Eqn. (|T^) that specifies the 
dynamics of one physical system has its J a depending on the fields of another, 
then we say that the second system interacts on the first. However, we must 
exercise some care in formulating this idea precisely. For instance, Ja is 
not gauge-invariant (and, indeed, can, by a suitable gauge transformation, 
be made to vanish); and furthermore it is unclear what "depends on" is to 
mean for fields on a bundle space. 

Consider again Eqn. (^. Fix a reduction of that system, so we have a 
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quotient bundle b ^ b ^ M, together with a field k/^a' on 6, satisfying 
conditions i) and ii) above. Any cross-section, of b gives rise to a new 
bundle, 6; and this 0, together with a cross-section of 6, specifies a cross- 
section of the original bundle b. Furthermore, the dynamical term k^^a' 
on b splits into corresponding dynamical terms on b and b. We say that the 
physical system represented by cross-sections of 6 interacts on the physical 
system represented by cross-sections of 6 provided there is no similar way 
to split the entire partial differential equation (|^) on b. Thus, in more detail, 
the 6-system interacts on the 6-system provided there exist no fields k^'^a, 
on b such that 

(VO/^i'". = /ii^^A^a, J A = l^A^JA, (27) 

up to gauge. This is to be compared with Eqn. (|25|). We are merely strength- 
ening condition ii) of the definition of a reduction to require an entire partial 
differential equation on b whose puUback is a linear combination of the given 
partial differential equation on b. In order words, if we can write the equa- 
tions on b in such a way that "no 6-fields are involved", then b does not 
interact on b. 

As an example, consider an electromagnetic field and uncharged perfect 
fiuid, in the presence of a background metric and derivative operator. Then 
the electromagnetic field does not interact on the perfect fiuid: We can find 
a reduction of this system in which the electromagnetic field is carried in 
b, the perfect fiuid in b; and a system of equations on b (the perfect-fiuid 
equations, described by k^"^si and j^) not involving the electromagnetic field. 
Similarly, the perfect fiuid does not interact on the electromagnetic field. 
However, if this is a charged perfect fiuid, then each of these physical systems 
interacts on the other (the perfect fiuid on the electromagnetic field through 
the charge- current term in Maxwell's equations; the electromagnetic field on 
the perfect fiuid through the Lorentz- force term in the fiuid equations). As 
a second example, note that the derivative operator interacts on virtually 
every physical system (through the use of this derivative operator in writing 
out field equations); and virtually every physical system interacts on the 
derivative operator (through Einstein's equation). 

We remark that the definition is manifestly gauge-invariant. Note that, 
as the definition is formulated, "interacts on" is not defined at all unless we 
have an appropriate reduction. Thus, for example, we are not permitted even 
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to ask whether the metric interacts on the electromagnetic field. (Perhaps it 
would be more natural to extend the definition so that background fields for 
a given physical system automatically interact on that system.) Finally, we 
remark that "interact on" need not be reciprocal: It is possible for system A 
to interact on B, but not B on A. It is not difficult to construct an example, 
e.g., with A a Klein-Gordon system and B a perfect fiuid. Introduce an 
additional term, involving the Klein-Gordon fields, on the right side of the 
equation, of fiuid particle-number conservation. Conservation of fiuid 
stress-energy is not thereby disturbed. However, I am not aware of any such 
examples arising naturally in physics. Should this observation be elevated to 
a general principle? 

We are, in a real sense, finished at this point. The world is described — 
once and for all — by the combined bundle B, with its cross-sections subject to 
(|T¥). In particular, whatever interactions there are in the world have already 
been included in the term of this equation. However, it is traditional 
to think of interactions in a somewhat different way — to think of them as 
capable of being "turned on and off' by some external agency (presumably, 
us). Consider, for example, the electromagnetic-charged fluid system. The 
flelds are Fab, n, p, and m"; and there appears, on the right side of the 
Maxwell equation (p9D, a term enu^ (charge- current), and, on the right side of 
the fluid stress-energy conservation equation (|iO| ) a term enF"'mU"^ (Lorentz 
force). Here, e is some flxed number (charge per particle). In the traditional 
view for this particular system, we think of this system as arising, not full- 
blown in its flnal form, but rather in two distinct steps. First, introduce the 
system with "no interaction" (e = 0), and then "turn on the interaction" by 
adjusting e to its correct value. 

This traditional view may be expressed, in the present general framework, 
as follows. On the combined bundle B, there is to be specified a "basic ver- 
sion" of the dynamical equations, ([T^) — a version in which "all interactions 
that can be turned off have been." Thus, the "electromagnetic interaction" 
has been turned off in the basic version: The only physical system that 
interacts on the electromagnetic field is the derivative operator. The "gravi- 
tational interaction" has been turned off: No field interacts on the derivative 
operator. For more complicated interactions — e.g., those of contact forces 
between materials — it may not always be clear how this "turning off' is to 
be carried out. The derivative operator generally survives into the basic ver- 
sion to interact with most physical systems. Physically, this phenomenon is 
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a reflection of tlie equivalence principle. Mathematically, it arises because it 
is difficult to eliminate the derivative operator and still lift diffeomorphisms 
(i.e., "maintain covariance") as in Sect. 5. This role of the derivative op- 
erator is related to the fact that the behavior of Va under an infinitesimal 
diffeomorphism, given in Eqn. (|16D, involves the second derivative of the 
vector field whereas other physical fields involve only the first derivative. 
To see this, consider a physical field, such as the Fat for electromagnetism, 
having s = 1 in ([T6|). Consider that portion of the combined equation, ([T^), 
referring to the electromagnetic field, and apply to it the infinitesimal dif- 
feomorphism generated by Then the term "(V$)m"" acquires a second 
derivative of So, invariance of ([I^) under infinitesimal diffeomorphisms 
can be maintained only if a second derivative of appears elsewhere in 
the electromagnetic portion of the equation, i.e., in the term J a- That is, 
some physical field must have, for its transformation behavior (|TBp under 
infinitesimal diffeomorphisms, s = 2. That field is the derivative operator. 

In addition to this basic version of ([l^), there is to be given some "inter- 
action fields", i.e., some fields 6Ja on B. We are free to add, to the left side 
of Eqn. ([I^) , any linear combination, with constant coefficients (the coupling 
constants), of these 6Ja- This is the step of turning on the interactions. One 
such 6Ja, for instance, is that which inserts the fluid charge-current nw" into 
(p9D, and the Lorentz force nF^mW™ into (0). Another inserts the trace- 
reversed combined stress-energies of all physical systems into0 Einstein's 
equation, (^). It is not entirely clear how much of this traditional view is 
psychological and how much physical. One way to argue that it has physical 
content might be to produce a general construction, given only the full com- 
bined Eqn. (0), that yields the basic version of this equation, as well as the 
appropriate 6Ja- 

These SJa do not span, at each point of the bundle B, the space of vectors 
of type "o"a" • That is, there are algebraic restrictions on the allowed interac- 
tions. These restrictions appear to be an essential part of the physical content 
of the systems under considerations. Disallowed, for example, are interac- 
tions on the electromagnetic field that represent a magnetic charge- current; 
on a perfect fluid that provide a source for particle-number conservation; and 

^''Note that these individual stress-energies cannot be inserted one at a time, with 
separate 6Ja, for to do so would violate integrability of the constraints of the combined 
system. 
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on the metric in p7|). See Appendix A for further examples. 

The basic version of Eqn. ([l^ on B must certainly be a viable system 
of equations, and so in particular its constraints must be complete[^ and 
integrable. What happens to viability of this system on turning on some 
interaction 5Ja? Completeness of the constraints will not change, for this 
involves only the dynamical term, that with coefficient Ka^^' in (0)- But 
integrability could be destroyed by turning on such an interaction. Indeed, 
the necessary and sufficient condition that integrability be retained for fixed 
constraint c"^" of Eqn. ( |14|) under addition to J a a term 6Ja is that, for some 

and Ta^ with Ta"" = 0, we have 

V„(c^'"(5J^) = (t^^^^Ja 

+r^(irA'"a + \iJA + SJA)iV7T)^n + r^(V7^),^ (28) 



where a a is the tensor that appears in ([TTD for the basic version of Eqn. 



(14). The proof of this assertion is straightforward: Demand that the inte- 
grability condition, (pT]), be satisfied for both Ka^o, J a and Ka^^w Ja + ^Ja- 
Thus, any viable candidate 5Ja for an interaction that can be "turned on" 
must satisfy the condition of (|2^). 

Eqn. (128|) is apparently a rather severe restriction on the interactions 
allowed in nature. The main reason for this is its nonlinearity: Given 5Ja 
and 5J'ai for each of which there exist r's satisfying (pSf), we have no guaran- 
tee that there exist r's that work for their sum. Here is a physical example 
of this behavior. Let 5Ja insert a perfect-fiuid stress-energy into Einstein's 
equation, and 5J'a insert an electromagnetic interaction into the perfect-fiuid 
equations. Then each of these modifications of (|n|) preserves integrabil- 
ity of the gravitational constraint, but their sum, 5Ja + 5J'a, does not. To 
achieve integrability in this example, it is necessary to add to this sum a 
term that also inserts the electromagnetic stress-energy into Einstein's equa- 
tion. Indeed, it is not even obvious, from (123), that whenever 5Ja preserves 



integrability in (|1^), then so does 2 6Ja^- 

So, finding a collection of interaction expressions, 5 J a-, that can be turned 
on in any linear combination, preserving all the while integrability of the con- 



straints of (14) does not appear to be easy. So, how does nature accomplish 



^^That is, complete in the sense appropriate to the difFeomorphism freedom inherent in 
dU) , namely that the left side of Eqn. (^ exceeds the right side by four. 
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this? Is there, for example, some simple, general criterion that can be applied 
to the 5Ja's to guarantee integrability? 



Appendix A — Examples 

Electromagnetism 

The field is an antisymmetric tensor field. Fab = F[ab], on M, with back- 
ground metric Qab- The equations are 

V'Fab = 0, (29) 

V[aFbc] = 0. (30) 

Thus, the fibres of the bundle b are six-dimensional, a typical vertical 
vector being given by 50" = 6 Fab = ^F^ab]- A typical vector in the space of 
equations is given by = (s",s"*'^), where is a vector (the coefficient of 
(p9D ), and s"*'^ an antisymmetric third-rank tensor (the coefficient of (^OD). 



Thus, the equation-space is eight-dimensional. 

The general hyperbolization at a point is given by 

h\,Sr' = iSF'^mt"', -U''5F'"^^), (31) 

where t°- is an arbitrary timelike vector. The general constraint at a point is 
given by 

^An ^ ^^gan^ ye"'"'''), (32) 

where x and y are numbers. Thus, the constraints form a two-dimensional 
vector space, while the space of vectors of the form c^"'nn, for fixed n„, is 
also two-dimensional. The constraints are complete and integrable. 
The allowed interactions are 

SjA={ja, 0), (33) 

i.e., no magnetic charge- current is allowed. These will preserve integrability 
provided V^ja = 0. 

Klein-Gordon 



34 



The fields consist of a scalar field, ip, and a vector field, ipa, on M, with 
background metric gab- The equations are 

= V'a, (34) 

"^[aA] = 0, (35) 
V>. = 0. (36) 

Thus, the fibres of the bundle b are five-dimensional, a typical vertical 
vector being given by 50° = {6tp,6ilJa)- A typical vector in the space of 
equations is given by cr^ = (s'*, s'^'', s) (the respective coefficients of (|3^-(|36D), 
where s'^'' is antisymmetric. Thus, the equation-space is eleven-dimensional. 

The general hyperbolization at a point is given by 

h^a'Sr' = i-W'S^, -tHtlj'\ ^t-6A), (37) 

where t"' is any timelike vector (say, future-directed), and w"' any vector not 
past-directed timelike or null. In order that this hyperbolization be causal, 
we must require in addition that w"" be future-directed timelike or null. The 
general constraint at a point is given by 

where x"" and y"'^'^ are arbitrary antisymmetric tensors. Thus, the constraints 
form a ten-dimensional vector space, while the space of vectors of the form 
c'^^rin, for fixed n„, is six-dimensional. The constraints are complete and 
integrable. It turns out that there are two separate physical systems (in the 
sense of Sect. 6) here. The first involves the field if) alone: The equation 
is (0), the hyperbolizations those generated by w"" in Eqn. (|37|) , and the 
constraints those generated by x""' in Eqn. (^). There are no background 
fields for this physical system. The other involves the field ipa'- The equations 
are (^), (0), the hyperbolizations those generated by in Eqn. ( 
the constraints those generated by y"*" in Eqn. (|38D. The background field 
is the metric. 

The situation regarding allowed interactions is not entirely clear (reflect- 
ing, perhaps, a certain lack of physical context for this system). Certainly, 
one allowed interaction is 

jA = (0, 0, -m'lP), (39) 
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where m is a number. This results in the massive Klein-Gordon system. 
(Note that the "tachyon equation" — the result of letting m be imaginary 
in ( PPP — admits a hyperbolization.) It seems likely that all allowed inter- 
actions have zeros in the first two entries of (0), for otherwise it is difficult 
to preserve integrability. (Passage to a charged Klein-Gordon system is not 
turning on an interaction in our sense, for it requires an entirely new bundle 
b. This example is discussed later in this Appendix.) Whether there are 
allowed other interactions of the form (|39D, but with different third entries 
on the right (for all of which, incidentally, all the constraints are integrable), 
is unclear. 



Perfect Fluid 

The fields consist of two scalar fields, n and p, and a unit timelike vector 
field, m", on M, with background metric Qab- The equations are 

{p + p)u"'VmU^ + {g^' + u^u'')VbP = 0, (40) 

«"V^P+(p + p)V^m"^ = 0, (41) 

u'^V^n + nVmU"" = 0, (42) 

where p{n,p) is some fixed function (the function of state). The first two 
equations are the components of conservation of T"* = (p + p)u"'u^ + pg"''', 
the third conservation of A^*^ = nu"'. 

Thus, the fibres of the bundle b are five-dimensional, a typical vertical 
vector being given by 50" = {6n,6p,6u'^), with Su" (because of unit-ness 
of m") orthogonal to -u*^. A typical vector in the space of equations is given 
by (j^ = {sa,s,s) (respective coefficients of (^OD-(^)), with Sa orthogonal 
to (reflecting that the left side of Eqn. (|^) is). The equation-space is 
five-dimensional. 

This physical system admits no hyperbolization unless 

P + J>>0. (|P + ^|P)>0. (43) 

op p + pan 

Physically, this is the requirement that inertial mass and sound speed both be 
positive. So, the fibres of the bundle b must be suitably restricted to achieve 
everywhere. The most general hyperbolization at a point is then given 
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by 



dp dp dp dp 

xiiip + p)^^+n-)Su., -Sp, -5p) 

+y(5n-^5p)(0, -1), (44) 

p+p p+p 



where x and y are any numbers with xy > 0, and where we have set Sp 



[^)6n + {^)6p. These hyperbohzations are all causal provided 



* + ^*<l (45) 

dp p + pdn 

(i.e., physically, provided the sound-speed does not exceed light-speed). If 
(^SD fails, then none are causal. There are no constraints (as follows, e.g., 
from existence of a hyperbolization and equality of the dimension of the space 
of fields and that of equations.) 

The allowed interactions (e.g., electromagnetic, contact-force, etc.) are 
given by 

jA = ija, J, 0). (46) 

That is, arbitrary sources are allowed in the equation of stress-energy conser- 
vation, but none is allowed in the equation of particle-number conservation. 



Gravitation 

The fields consist of a symmetric, Lorentz-signature metric, Qab, together 
with a (torsion-free) derivative operator, Va, on M. The equations are 

"^aQbc = 0, (47) 

Rab{c^9d)m = 0, (48) 
Rmiab)"" = 0. (49) 

Here, Rabc'^, with symmetries Rabc'^ = R[ab]c'^ and R[abcf = 0, is the curva- 
ture tensor (i.e., the "derivative of the derivative operator"), defined by the 
condition that 

V[,Vfe]ec = ^Rabc%, (50) 

for every covector field on M. 
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Thus, the fibres of the bundle b are fifty- dimensional, a typical verti- 
cal vector being given by 6(j)°' = {Sgab,ST°-bc), where Sgab = 5g(ab) (ten di- 
mensions), and ^r^bc = '^r"(6c) (forty dimensions). The latter represents a 
first-order change in the derivative operator, whose effect is to replace V a^b 
by V aib + ^^^abim- A typical vector in the space of equations is = 
^^abc^^abcd^^ab^ (the respective coefficients of (|43)-(||)), where s"^^ = s"^''^). 



gabcd _ g[ab](cd)^ ^ab _ g(ab) ^ Thus, the cquatiou spacc has dimension one 
hundred ten (= 40 + 60 + 10). 

This system consists of two separate physical systems (in the sense of 
Sect. 6). For the first, the field is the metric gab, and the equation (^Tj), with 
no background field. Its most general hyperbolization at a point is given by 

h^o^'Sr' = x''''"'^6g^n, (51) 

where x'^^'^™"' has symmetries x"'"^™'" = x"-^"^"'^^'"^\ and is such that UaX"''"^"^^ 
is positive-definite in the index pairs "6, c" and ^^m,n" , for some n^. The 
general constraint for this physical system at a point is given by 

where x""'"^ = x^"'°''^^'"^\ Thus, the vector space of constraints has dimension 
sixty, while the space of vectors of the form c^"'nn, for fixed n„, has dimension 
thirty. These constraints are complete and, by virtue of Eqn. (^8]), integrable. 



For the other physical system, the field is the derivative operator, the equa- 
tions are (|^), (^), and the background field is the space-time metric gab- 
This physical system has no hyperbolization, because of the diffeomorphism 
freedom. (But it does have a "hyperbolization up to diffeomorphisms" . See 
Eqn. (ESf)). The most general constraint for this system is given by 



^An ^ (^^nabcd ^ 2x^'' g^^^^ g^> - g^'^ x^"" g^^"" , x^^g^^"" - ^^''''x"), (53) 

where x""*'^'^ = xt"'"^]*-'^'^''. Thus, the vector space of constraints has dimension 
fifty-four, while the space of vectors of the form c^^'rin, for fixed n„, has 
dimension thirty-four. These constraints are integrable, but not complete. 
(This feature, again, is related to diffeomorphism-invariance. See Sect. 5.) 
The most general allowed interaction, apparently, is 

6jA = {0, 0, Tat-^T'^mgab), (54) 
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where Tab = Ti^ab) (the stress-energy of matter). In order that this inter- 
action preserve integrabihty of the constraints, we must have conservation, 
V^Tab = 0. This is a very severe restriction on the fields contributing to Tab, 
and the interactions between those fields. 

Spin-s Systems 

The field is a totally symmetric, 2s-rank spinor, ip^'"^ = ip'^^"'^\ on M, 
with background metric gab- Here, s = |,l,|,2,---. The equation is 

V^'^^^^-^ = 0. (55) 

Thus, the fibres of the bundle h are (4s -|- 2)-dimensional. (This, and all 
subsequent, reference to dimension means real dimension.) A typical vertical 
vector is given by 5(j)^ = S^"^'"^ = Sip^^"'^\ A typical vector in the space of 
equations is givenQ by = sa'b---d = sa'(b-d)- Thus, the equation-space 
is 8s-dimensional. 

The most general hyperbolization at a point is given by 

h^a'Sr' = tB'-D'B-DSijA''''-''', (56) 

where tB'-.-D'B - D = i(B'-D'){B - D), a Hermitian quadratic form on symmetric, 
(2s — l)-rank spinors, is positive-definite. The general constraint is given by 

^An ^ ^N'A'^N(B^C-D)^ (57) 

where x^'"^ is an arbitrary symmetric, (2s — 2)-rank spinor. Thus, the con- 
straints form a (4s — 2)-dimensional vector space, while the space of vectors 
of the form c"^"n„, for fixed n„, also has dimension (4s — 2). Note that, for 
s = 1/2, there are no constraints. These constraints are always complete, but 
they are integrable if and only if either s < 1 , or the metric gab is conformally 
flat. (This is the famous "inconsistency of the higher-spin equations".) Ex- 
cept for the cases s = 1 (electromagnetism), and s = 2 (linearized gravity), 
it isn't clear what are the allowed interactions. For charged spin-s fields, see 
later in this Appendix. 

Elastic Solid 

^^The notation is a little awkward here. The index on the left is in the space of equations; 
those on the right, in spinor space. 
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The fields consist of a scalar field p, a unit timelike vector field m", and a 
symmetric tensor field hat of signature (0, +, +, +) satisfying habu'' = 0, with 
background metric gab- The equations are 

u'^V^hab + 2{ViaUnhb)^ = 0, (58) 

pu^VmU'' + qW^r""' = 0, (59) 

u'^VmP + pVmU'^ + r'""V™n„ = 0, (60) 

where we have set qab = gab + UaUb, the "spatial metric" . Here, r"^ is some 
fixed algebraic function of hab, m", and gab, satisfying r^^-Ub = and r"* = 
^(ab) rjnj^^ physical meaning of these equations is the following. The field p is 
the mass density, and u"" the material four-velocity. The field hab represents 
a sort of "natural spatial geometry" for the material. Thus, we interpret the 
combination hab — Qab (natural geometry minus actual geometry) as the strain 
on the material; and Eqn. (^) (which is just Cuhab = 0) as requiring that the 
material carry along with it its natural geometry. The field r'^* represents 
the stress of the material. Thus, we interpret T°'^{hab) as the stress-strain 
relation^, and the combination pu°'u^ + r"^ as the material stress-energy. 
Eqns. (0), (^) are precisely conservation of this stress-energy. 

This system admits a hyperbolization if and only if p > 0; and in addition 
the tensor r"'^'^'^ = Ot^^ / dhf,c is symmetric under interchange of the index pairs 
a,b and c,d and positive-definite in these pairs. In this case, the most general 
hyperbolization is given by 

h\,Sr' = x{25p-h'^"'{pq^^ + T^m)Sh,,]{-h'^^{pq\ + T\), 0, 2) 

+y{T-'-''5Kd, 2h-^5um, 0), (61) 

where x and y are numbers with > 0, and h""^ = h^"'''^ is defined by 
Ubh"-^ = and h"'"^hbm = q°'b- These hyperbolizations are all causal if and 
only if 

/^rnnr"'""" < ^pg'^^ (62) 

which means, physically, that no acoustic wave-speed exceed the speed of 
light. This system has no constraints. If we generalize this system to allow 

^^Note that the stress-strain relation need not be Hnear. It would perhaps be natural 
to require that r"'' vanish when hab = <lab: but that requirement is not needed for what 
follows. 
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the stress r'^^ to depend, not only on the natural geometry hat, but also on 



the mass density p, then the Eqns. (58)- (60) never admit a hyperbolization. 



It seems peculiar that such a physically benign generalization would preclude 
a hyperbolization. What is going on here? 

Presumably, the most general allowed interaction is given by 

JA=(0, Ja, j). (63) 

That is, interactions are allowed that exchange energy-momentum with the 
environment, but not that modify Lie transport of the natural geometry. 

Special Relativity 

The fields consist of a symmetric, Lorentz-signature metric Qab, together 
with a derivative operator Va, on M. The equations are 

Vagbc = 0, (64) 

Rab/ = 0, (65) 

where Rabc^ is the curvature tensor, given by (|50|) . 

The bundle space b is this case is identical to that for general relativity, 
and so in particular the fibres have dimension fifty. But now the equations 
are different. A typical vector in the space of equations is = [s"-'"^, s°''"^d) 
(respective coefficients of (0), (|65D), where s"'"^ = s^-'^^d = s'^^'l'^rf, and 

g[abc]^ _ Q TJius^ the equation-space has dimension one hundred twenty 
(= 40 + 80). 

This system consists of two separate physical systems (in the sense of 
Sect. 6). The first, with field the metric, is identical to the similar system 
for general relativity. Thus, the hyperbolizations are given by (^TJ), the con- 
straints by (|^). For the second system, the field is the derivative operator, 
and the equation ([651) . This system has no hyperbolizations, because of the 
diffeomorphism freedom. The most general constraint for this system is given 
by 

with x'^"''"^d = x'""^''^^ and xt"'"^'^'^ = 0. Thus, the vector space of constraints 
has dimension sixty, while the space of vectors of the form c^^n^, for fixed 
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has dimension forty-eight^ The constraints are integrable, but not complete. 

Apparently, no interactions whatever are permitted in Eqns. (|64D , (|65|) . 
Note that passage to general relativity is not "turning on an interaction", 
because this is a change in the dynamical part of Eqn. (|65D. 

Dust 

The fields consist of a scalar field, p, and a unit timelike vector field, m", 
on M, with background metric gab- The equations are 

u'^VmU^ = 0, (67) 

M^V^^p + pVmU"" = 0. (68) 

Thus, the fibres of the bundle b are four-dimensional, a typical vertical 
vector being given by 50° = {6u"',6p), with UaSu"" = 0. A typical vector in 
the space of equations is = {sa,s) (respective coefficients of (167|) , (|68|)), 
with SaU"' = 0. The equation-space is four-dimensional. 

Remarkably enough, this system admits no hyperbolization. To see this, 
note that the most general candidate for a hyperbolization at a point is 

/i^^,50°' = {xabSu^ + Xa5p, yaSu" + y5p), (69) 

for some tensors Xab, Xa, ya, and y. Combining this with the kA^a' from 
Eqns. (0), (0), we obtain 

/i^„/50°'fcA™/3'50^' = u"'[xabSu''6u'' + Xa^u^Sp + VaSpSu" + y5 p5 p] 

+p5u^[yju'' + y5p]. (70) 



We see that this is symmetric in 50° and 50^ if and only if ya-, and Xa 
all vanish, and Xab is symmetric. But these conditions preclude positive- 



definiteness of (|70|). This lack of a hyperbolization is discussed briefly at the 
end of Sect. 3. There are no constraints. 

Any 5jA is, apparently, allowed as an interaction. 



-"^^Here is a mystery. For this system, the left side of Eqn. (g) exceeds the right side by 
eight. But we might have expected, from the fact that Eqns. ([64[), ( |65| ) have an initial- 
value formulation up to diffeomorphisms, an excess of four. What is the explanation for 
this discrepency? 
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Charged Fields 

Fix on M a smootli, antisymmetric tensor field Fab having vanishing curl 
(V[a-^6c] = 0). We assume that M is simply connected, and that F has been 
so chosen that its integral over every compact 2-surface in M vanishe^^ 

We wish to introduce the notion of tensors on M of charge e (some fixed 
number) and index type t°'"'b--- (some fixed arrangement of contravariant and 
covariant M-indices). To this end, fix a reference point Xo of M, and consider 
any point x of M. Consider the collection of all pairs, {Aa,t"'"'b-.-), where Aa 
is a vector field on M with V[aAb] = Fab (existence of which is guaranteed 
by the assumptions above), and t°'"'b- - is a complex tensor at x. Two such 
pairs, {Aa,t°'"'b-. ), and {A'^,t"^"'b-. ), are taken as equivalent if 

f'^ -b- = exp[te J {A'^ - A^)rfs™] r-b..., (71) 

where the integral on the right is over any curve from Xq to x. (Independence 
of the choice of curve is guaranteed by the conditions, above, on Aa and 
A'^.) This is an equivalence relation. The equivalence classes, called charge-e 
tensors at x, form a vector space of (complex) dimension 4*, where s is the 
total number of indices of f^ ' h-- A charge-e tensor field is an assignment, to 
each point x G M, of a charge-e tensor at x. Thus, a charge-e tensor field can 
be represented by a pair {Aa,t"-"'b-. ) of fields on M, where pairs {Aa,t"'"'h-) 
and {A'^.t'"-'" b-) are identified provided Eqn. ([TTl) holds for all x. 

The discussion above was predicated on the choice of a fixed reference 
point, Xo G M. Let us now change to a new reference point, Xo G M. 
Fix a smooth curve 7 from Xo to Xo- Then we identify a charge-e tensor 
at X, defined via reference point Xo, with charge-e tensor at x, defined via 
reference point Xo, provided these have respective representatives, (Aa, t°'"'b--) 
and (Aq, t"'"b...), with 1"-"^... = exp[ze A^iis™'] t°-"'b---- Note that this is 
independent of representative. But it does depend on choice of the curve 7: 
A change in 7 changes this identification by an overall phase (the same for 
each point x). Thus, charged tensor fields (independent of reference point) 
make sense only up to an overall constant phase. 

^^Vanishing of the 2-surface integrals of F means physically that all wormholes manifest 
zero net magnetic charge. In fact, we could allow here integral multiples of a certain 
fundamental magnetic charge, but, for simplicity, do not. The assumption of simple- 
connectedness of M avoids our having to consider Aharonov-Bohm effects. Again, we 
could relax this assumption, at the cost of somewhat complicating the discussion below. 
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Now let there be given on M a space-time metric gab, with corresponding 
derivative operator Va- We extend the action of this derivative operator to 
charged fields as follows. Given a charge-e tensor field, choose any repre- 
sentative {Aa,t°-"' b-. ), of it, and let its derivative be the charge-e tensor field 
with representative {Aa, (Vc — ieAc)t°'"'b-), noting that this is independent 
of the original choice of representative. 

The charged Klein-Gordon system consists of charge-e fields ip and ipa, 
with background fields the metric gab and electromagnetic field Fab- The 
equations are the same as (p4|)-(^), where, of course, Va now is the derivative 
operator on charged fields, except that there appears on the right of Eqn. (p5|) 
the term ieipFab- (Failure to include this term would destroy integrability of 
the constraints.) Similarly, the charged spin-s system consists of a charge-e 
spinor field, ip^'"^, satisfying Eqn. ([551) . 

In order to fit these systems into the present framework, we must fix a 
reference point Xo, for otherwise the freedom to multiply fields by an overall 
constant phase makes it impossible to find a bundle to house the charged 
fields. The discussion of hyperbolizations and constraints then goes through 
in a manner identical to that of the uncharged case, with one exception. The 
constraints for the equation for charged spin-s fields, for s > 1 and Fab 7^ 0, 
fail to be integrable. One further complication arises. There is no natural 
way to lift diffeomorphisms from M to the bundle spaces b associated with 
any of these charged fields (although there does, of course, exist a lift "up to 
overall constant phase"). The reason for this is our having fixed a reference 
point, Xo, in order to introduce b. 

Kinetic Theory 

Let M, gab be a time-oriented space-time. Fix a nonnegative number m, 
and denote by F the seven- dimensional manifold of all pairs (x, pa), where 
X G M, and Pa is a future-directed vector at x with paP"' = —m^. Thus, F 
is a fibre bundle over M. Denote by F^. the fibre over x G M, so each F^. is 
a three-manifold. The field for kinetic theory is a nonnegative function / on 
the manifold F; and the equation is 

p'^V./ = C(/). (72) 

Here, C is, for each x G M, a mapping from nonnegative functions h on F^; 
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to functions on T^, satisfying 

J C{h)p'' = 0, y C(/i)py = 0, (73) 

for every h, where the integrals are over using the natural volume element 
on this mass shell. Eqn ( |72D is to hold for each (x, Pa)^ T, with C evaluated 
on the restriction of / to Fa;. The physical interpretation of these equations 
is the following. The nonnegative function / is the distribution function (of 
particle position-momentum) for mass-m particles, and Eqn. ( [72D is Boltz- 
mann's equation. The C in ( [7^ ) is the collision function; and ( [73| ) is local 
conservation, in collisions, of particle-number and energy-momentum. 

This system does not, strictly speaking, fall within the famework of Sect. 
2, for the space of allowed "field values" at each point of M (nonnegative 
functions on T^) is infinite-dimensional. But, if we agree to ignore this one 
defect, there results a nice example of our framework. A typical vertical 
vector is given by 50" = 5/, a function on Fa-. A typical vector in the space 
of equations is o"^ = s, a function on F^,. The most general hyperbolization 
of (|7^ is given by 

= /i 5/, (74) 

where /z is any positive function on F^,. There are no constraints. Presum- 
ably, there is allowed as an interaction in (|72]) any function 5j on F such 
that, for every x G M, J{Sj)p°' = 0, where this integral is over F^.. That 
is, interactions may not violate local particle-number conservation, but are 
otherwise arbitrary. 

Lagrangian Systems 

Fix a fibre bundle b ^ M. By a (first-order) Lagrangian on the bundle 
6, we mean a smooth function L as follows: L is a function on pairs, (k, 
Ca"), where k is a point of the bundle space 6, and Ca" is a tensor at k 
satisfying Ca"(V7r)a^ = ba\ and, for each such pair, L^nXa"") is a density^ 
in M at the point 7r(K). Such a Lagrangian gives rise to a system of partial 
differential equations on cross-sections of the bundle b. In order to write 
these equations explicitly, it is convenient to introduce a derivative operator 
Vq, on mixed fields on b, such that the derivative of every vertical vector field 
is vertical, and VQ(V7r)^" = 0. Then, e.g., the operator "derivative along the 

^"^That is, an antisymmetric, fourth-rank Af -tensor, whose indices we shall suppress. 
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cross-section 0" is given by (V0)a"Va. (See the discussion just preceeding 
Eqn. ([l^)-) Written in terms of this operator, Lagrange's equation becomes^ 

57(V0),^^V^((V0)/) + (V0)^^V^(-— ) - V,,L = 0. (75) 



The coefficients in ([75| ) are to be evaluated at Ca" = (V0)a". This equa- 
tion is of course independent of the choice of derivative operator. It is also 
unchanged under adding to L any function of the form (Vof ")Ca" (a "total 
divergence"), where v"" is any M-density field on h. Note that the spaces of 
equations and unknowns for Eqn. ( |75D have the same dimension, namely, 
that of the fibres of the bundle h. 

In order to cast Eqn. (^) into first-order form, we introduce auxiliary 
field Ca", subject to Ca'^(V7r)Q* = 5a ] and we supplement ( [75| ) with the 
additional equations 

(V0)." = Ca", (76) 
C[/V|^|Ca]" = 0. (77) 
The system ([7BD, ([77|), ( |75D is closely analogous to the Klein-Gordon system, 

m-m- 

The fibres of the bundle h appropriate to this system have dimension 5n, 
where n is the dimension of the fibres of h. A typical vertical vector in h is 
given by = (50" , ^Ca" )■ Here, 50" , a vertical vector in 6, represents 
an infinitesimal change in the value of the cross-section of 6, while dd^a" 
represents an infinitesimal change in the tensor (^a". A typical vector in the 
space of equations is = (s^^/, s"*o', s" ), with|^ s'^^a' = st^'^'a/ (respective 
coefficients of (|76D , ([f7|), (|75D) Thus, the equation-space has dimension lln 
(= An + Qn + n). 

Set 

the coefficient of the first term in (|75|) . The system ([751)- ([77D admits a hyper- 
bolization at a point if and only if the tensor S°'^a'i3'V°' is Lorentz-signature 

^^Note that the Greek index of "d/dCm" " acquires a prime, as a consequence of the 
fact that Ca" (V7r)„'' = (5a^ 

^^Note that the difference between the two sides of (^6|), as well as the left side of (|T 
is automatically vertical. Thus, a prime must be attached to the Greek subscripts of s° 
and s^'^a'- 
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for every nonzero f"', and, in addition, there exists a tangent vector t° and 
a covector Ua that are both timehke for every one of these Lorentz metrics. 
Note that this is a rather severe condition on S'^'^^'p'- When it is satisfied, 
the most general hyperbohzation is given by 

= {-w\.p>, -t^'^S'^^^^p^dCJ', ^t'^^C™"'), (79) 

where t"- is such a common timehke vector, and w"'a'/3' = w°'{a'i3') has the 
property that {t^nf))naW°'a'i3' is positive-definite for some such common time- 
hke covector n^. (Such a w^^a'/^' always exists, e.g., that given by t^'Ga'p' 
with Ga'p' positive-definite.) Given a space-time metric gab-, this system is 
causal if and only if every lying in one half of the light cone of Qah will 
serve as a "common timelike covector" above. Eqn. (^) should be compared 
with Eqn. (p7|), giving the Klein-Gordon hyperbolizations. The most general 
constraint for the system ([75D-(|77|) is given by 



c^" = (x"V, l/"'^^', 0), (80) 

where x'^°'a' = x^'^'^^a' and y"''^^a' = y^^"'^'^a'- Thus, the constraints form a 
vector space of dimension lOn, while the space of vectors of the form c^"'nn, 
for fixed n„, has dimension 6n. Eqn. (|80|) should be compared with Eqn. 
(|38D, giving the Klein-Gordon constraints. These constraints are complete 
and integrable. 

This system only allows those interactions that preserve its Lagrangian 
character, i.e., that result from some change in the Lagrangian. This change 
must be so chosen to leave invariant the tensor S^-^a'/s' of Eqn. ( |78D (this being 
the coefficient of the dynamical term in (^) ) . The most general such change 
is given by 6L = W^aCrn" + W , where W^^ and W are arbitrary smooth 
fields on the bundle space h. That is, these W^s are allowed to depend on 
the K- variables, but not on the Ca"- This change in the Lagrangian L results 
in a term 

SjA = (0, 0, 2{y^^Wp^^)CJ - Wc^'W) (81) 

in (0). Thus, the allowed interactions on a Lagrangian system are very special 
indeed. In particular, the 5jA can be at most linear in the field Ca"- 

One example to which the discussion above can be applied is the Klein- 
Gordon system, with Lagrangian L = ^(Va'ip)(V"'4')- But there exists an al- 
ternative Lagrangian formulation for this system, starting from L = —\ipa'>P°'+ 
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ijj^'Va'^- This alternative Lagrangian involves the full set of Klein-Gordon 
fields, i^a (as opposed to just ip), and yields equations that are automat- 
ically first-order (as opposed to second-order). It turns out that such an 
alternative Lagrangian formulation is available quite generally. We sketch 
below how this comes about. 

Recall that the points of the bundle b are denoted k, and that a La- 
grangian on b consists of a certain function L{k, Ca")- Fix such a Lagrangian. 
Consider now the bundle b, whose points are pairs (k, Ca")- Let us now in- 
troduce, on this new bundle b, the following Lagrangian: 

L = {^)iy<P)a--L. (82) 

The Lagrange equations arising from (|^) are precisely (|75|) and (|76D , except 
that in ( ffSD we must replace "(V0)a"" everywhere by "Ca""- Thus, we ob- 
tain from the Lagrangian (|8^) a first-order system from the start. But this 
system — despite the fact that its spaces of unknowns and equations have 
the same dimension — admits no hyperbolization. In fact, this system has 
constraints — essentially those given by the x'^°'a' in Eqn. (PU[). These con- 
straints are not integrable, but their integrability conditions are (|77D , which 
are linear in first derivatives of the fields (as opposed to quadratic, which 
is what happens in the general case, (0)). So, we may include these inte- 
grability conditions with the other equations of our system. The result is 
Eqns. (|7|)-(0), a system that admits a hyperbolization and has complete, 
integrable constraints. That is, the result is a system with an initial-value 
formulation. 

Higher-Order Systems 

It is conceivable that some physical phenomena may be described by 
higher-order systems of partial differential equations (e.g., arising from a La- 
grangian of higher order) . We describe briefly the conversion of such systems 
to first-order form, and their resulting initial-value formulation. 

Consider a quasilinear, s^'^-order system of partial differential equations, 
which we may write in the form 

kA'^'-'^^a'Vm, ■ ■ ■ V^,r' + J A = 0, (83) 

where /ca*"^ ^^o' = kA^^^'"^"'' a' and ja are functions of point of M, the field 
0, and, at most, its first (s — 1) derivatives. To achieve this form, we have 
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introduced a connection in the bundle of which is a cross-section. (The 
coefficient k in (|83|) , but not j, is independent of this choice.) As an example, 
Lagrange's equation for a higher-order Lagrangian takes the form (^5]), with 
the index "A" replaced by with kp/"^'^'"'^'' a' symmetric in (3', a', and 

with s even. Let us now, in order to achieve a first-order system, introduce 
auxiliary fields, 0^" , (pab"' ■ 4>ai-as-i°' ; each symmetric in its Latin indices, 
and their corresponding equations, 

V[a0ai]a2-ai" = l^'aai-a"' ■ (85) 



Here, i = 1, ■ ■ ■ , (s — 1), and the /i on the right in (|85|) is a certain function 
of the lower-order 0's. Our ffist-order system consists of Eqn. (^) with 
the derivative-term replaced by "Vmi^ma - ms" '\ and Eqns. (|84D , (185|) . The 
constraints for this system are always complete, and, by virtue of the choice 
of /i in (^5]), integrable. 

When does this system admit a hyperbolization? No simple, general 
criterion is known, but the following remarks will at least suggest a possible 
line of attack. First note that the equations on 0ai-- ai° for i < {s — 1) 
(namely, (|85D for z < (s — 1), and (0)) always admit a hyperbolization (in 
a manner similar to that of the Klein-Gordon system, (|3^) and (^)). What 
remains is the field 0ai --as_i" , and its equations, ( ^5]) (for i = {s — 1)) and 
(|83D. Let t-^a' be any tensor such that t-^a'kA^^ (}' is symmetric in 
and let t"^"'"^-^ be any totally symmetric tensor. Consider the tensor given 
by 

p mi - ms_iana---ns_i _ ^ +A j.{mi - ms-i r, ni---na-i}a 
-Tq,' P' — -t a'l^ l^A 13' 

-it^„,t"^"^^-'"=-2A;A'"=-^"i-"^-i>^,, (86) 

where "{mi ■ ■ ■ms_ini ■ ■ -Us-iY^ means "add all 2(s — 1) terms that result 
from cyclic permutations of these indices, and then symmetrize the result 
over mi ■ ■ -m^-i and over rii ■ ■ -n^-i". Then, as is easily checked directly, 
this P has the properties 

p m\---ms-ian\---ns-i p {ni---ns-i)a(mi---ms-i) fQ'7\ 

^a' 13' — -t^{a' /3'), [oi) 

-Ta' 13' — t t a'l^A f3' ■ [oo) 
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It follows from Eqn. ( ]87| ) that the differential operator p^,"T-i- ms-i«"i-"s-i^, 
^a(pni--ns-i'^ IS automatically symmetric; and, from Eqn. (p8|), that this 
differential operator is a linear combination of the left sides of ( P^D and 
(|83D. Thus, we obtain in this way a symmetrization of our first-order system. 
This symmetrization is actually a hyperbolization provided we can choose 
t^a' and t°'^- -°-=-^ such that there exists at each point a covector for which 
the quadratic form naPa'"^^"'"^"~^°'"'^"'^"~^ fs' is positive-definite (i.e., is positive 
when applied to any nonzero 50mi-- ms_i° )• 

When can this positive-definiteness condition be achieved? For s = 1 (i.e., 
when (p3| ) is already first-order), this condition just repeats the definition of 
a hyperbolization of a first-order system. The complete solution for s = 2 is 
given in the discussion of Lagrangian systems in this Appendix. The cases 
s > 2 are considerably more difficult. Those for odd and even s appear to 
be rather different in character. Are there any simple, general conditions on 
j^^mi-nis ^ for s > 2, that imply existence of t'^a' and t"-^-"-"-^ yielding, by the 
construction above, a hyperbolization? Are there any other symmetrizations 
of the first-order system (p3D-(p^)? 

It is very likely (depending in part on how the derivative in (|8^) is struc- 
tured) that these higher-order fields would, under infinitesimal diffeomor- 
phisms, pick up in (|16D higher derivatives of the generating vector field 
Thus, as discussed in Sect. 5, such fields may usurp the diffeomorphism 
gauge group from general relativity. 

Appendix B — Existence and Uniqueness of 
Solutions of Symmetric, Hyperbolic Systems 

Fix a smooth, four-dimensional manifold M, and a finite-dimensional vector 
space F (tensors over which will be denoted, respectively, by Latin and Greek 

indices). We are interested in F- valued functions on M, M F. Consider, 
on such functions, the partial differential equation 

ka"'f3{x,(f)^x))Vm<P^ +Ux,r{x)) = 0, (89) 

where x G M, and ka"^/3 = k(^a'^/3) and ja are smooth functions on M x F. 
This is just Eqn. (0), modified as follows: i) we have written the bundle space 
as a product, and imposed a vector-space structure on the fibres (which can 
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always be done locally); ii) we have chosen the gauge so that the last index 
of k is vertical, and then, since now all Greek indices are vertical, have 
suppressed all primes; and iii) we have multiplied Eqn. through by a 
suitable hyperbolization, as in (|^). By initial data for (|89D, we mean a three- 
dimensional submanifold S of M, together with a smooth mapping S ^ F, 
such that, for every point x of S, the tensor n^, ka^i3ix,4>o{x)) is positive- 
definite, where normal to S at x. (This agrees with the definition of 
Sect. 3.) 

The fundamental theorem]^ on existence and uniqueness of solutions of 
symmetric, hyperbolic partial differential equations may now be stated as 
follows. 

Theorem. Let 5 ^ F be initial data for (H). Then: 

a. For some open neighborhood U of S, there exists a smooth solution, 
U^F, of (|H|), with 015 = 00- 

b. For U F and U' ^ F two such solutions (so 0|5 = 0'|5 = 0o), there 
exists a neighborhood ?7 C f/ fl f/' of S* in which = 0'. 

c. Any U G U n U' will serve for part b provided it can be covered by 
a smooth, one-parameter family, St, of three-sub manifolds of U such that i) 
one of the St is S itself, ii) each St coincides with S outside of a compact 
subset of S, and iii) for every t, both and 0', restricted to St, are initial data. 

Part a of this theorem is existence: It guarantees a solution of (p9|), 
satisfying the given initial conditions, in some neighborhood of the initial 
surface S. Part b is uniqueness: It guarantees that two solutions, each defined 
in some neighborhood of 5*, must coincide in some common subneighborhood. 
If. Part c strengthens part b by guaranteeing a certain minimum "size" for U : 
Part b holds for any U that can be covered by a one-parameter family, St, of 
surfaces that result from deforming S within compact regions, provided each 

St can serve as an initial-data surface^. Thus, given U F and U' F, 

^^F. Johns, Partial Differential Equations, Springer- Verlag (New York), 1982; P.D. Lax, 
Comm Pure and Appl Math 8, 615 (1955); K. O. Friedrichs, Comm Pure and Appl Math 
7, 345 (1954). 

^''This U lies within the domain of dependence of S, suitably defined. It is merely for 
convenience that we characterize U using deformations of S, rather than curves whose 
tangents are propagation directions, as defined in Sect. 3. 
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we can generate a U that works for part b by taking compact deformations 
of S within the intersection U D U', stopping as soon as positive-definiteness 
of rim ka"^ /3{x, (j){x)) OT Tim p{x , (j)' {x)) fails. It foUows in particular from 
part c that the solution at a point of this tj depends only on the data in a 
certain compact region of S. This observation is the basis of our discussion 
of "signal-propagation directions" in Sect. 3. We remark that there are also 
results that strengthen part (a), by guaranteeing a certain minimum "size" 
for its neighborhood f/, but these are more complicated and less useful than 
the strengthening of part b, as given by part c above. Of course, solutions of 
Eqn. (|89D can — and often do — evolve to become singular. 
A sketch of the proof of the theorem follows. 

Fix an alternating tensor, eabcd = ^[abcd\, on M (to facilitate integration 
and to allow us to take divergences of vector fields), and a positive-definite 
metric G^/? on F (to facilitate taking norms). We now derive an inequality, 
(|9TD, that we shall use three times in what follows. We first note that, for 
any fields 0", 0'", k^'^p = k^^J^p), h!j^p = fc'(„™/3), and on M, we have 

Vm[(ka"'f3+K'^p){r - 0")(0^ - 0'^)] 
= [V™(fc„™;3 + k'j^(3)W - - O 

-4(0°-r)[L-J'j 

+4(0" _ 0-)[fc„™^V„0^ + J, - A;L>V^0'^ -I], (90) 

where we have everywhere suppressed the variable x. To prove Eqn. (|90|) , 
expand the left side, and note that all terms cancel. Next, let St (t G [0,to]) 
be a smooth family of three-submanifolds of M, each of which coincides with 
S = So outside a compact subset of S, and on each of which Umika'^p + k^'^p) 
is positive-definite, where is the normal to St in the direction of increasing 
t. Denote by V the union of the St, so V has a boundary consisting of S 
and St^- Now multiply (|90|) by exp(— 2t/r), where r is a positive number, 
and integrate over V. Integrating the left side by parts, the resulting volume 
integral involves the expression ^(Vmt){ka™'/3 + k^"^/3). But, by construction, 
this tensor is positive-definite, and so this volume integral will, provided r is 
chosen sufficiently small, dominate the integral of the first term on the right. 
For the remaining three terms on the right in Eqn. (|90D , use the Schwarz 
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inequality. There results 

■J s 

+2ar| |(^ - V(0 + 0')e"*/"| | + 4(Tr| | (j - l')e-''^\\ 
+4ar 1 1 ( W0 + J)e-*/" 1 1 + 4ar 1 1 (^'V0' + J)e"*/" 1 1 , (91) 

where we have suppressed indices. The first term on the right in (^) is a 
surface term arising from the integration by parts (the surface term at St^ 
having been absorbed into the inequality). In (pT]), || || means "the square 
root of the integral of the square of the indicated field over (i.e., the 
L^-norm), and r and a are positive constants such that 

rV™(fc„™/3 + k'j-p) < {Vmm^'^p + :fcL™/3) > -G^p, (92) 

a 

everywhere in V . (Note that such constants exist, by compactness of V and 
positive-definiteness of the middle expression in (|92D.) The inequality (pT| ) is 
our final result. It asserts that two fields, and 0', are close to each other 
in V (left side of (pi])) provided that they are close on the initial surface 
S (first term on the right), that each approximately satisfies an equation of 
the form ( P5| ) (fourth and fifth terms on the right), and that their respective 
coefficients, (/c, j) and (fc , j ), in this equation are close (second and third 
terms on the right). It is in this derivation of the inequality (^Tj) that we make 
crucial use of the symmetric, positive-definite character of the coefficients 
kcJ^jS, and the geometrical conditions of part c of the theorem. 

We first prove uniqueness (parts b and c of the theorem). Let and 0' be 
two solutions, as in part b, and let U G U OU' and the family St be as in part 
c. Apply inequality (|1|), with fc^*^^ = koTaix, 0(x)), k'j'^p = kj^p{x, 0'(x)), 
3a = Ja(a;, 0(a;)), and = ja(x,0'(x)). Then the first term on the right 
vanishes (by initial conditions), and the last two terms on the right vanish 
(by (|89|)). But, for the remaining two terms, each of |fc(0) — fc(0')| ^ind 
|j(0) — j(0')| is bounded by a multiple (namely, the least upper bound of 
|(9/c/90| and |9j/50|, respectively) of |0 — 0'|. So, choosing r sufficiently 
small, the sum of these two remaining terms on the right of (|9l|) is less than 
the left side of this inequality. We thus conclude that ||(0 — 0')e~*/'^|| = 0, 
and so that = 0'. 

We next turn to existence. This is carried out in two steps. 
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Consider first the equation 



ka"'f3Vm^''+J^ = 0, (93) 

with fixed fields ka^pi^x), This is just (P^D, but with the coefficients 

evaluated on a fixed background field. Fix a family St of surfaces, and a 



region V , as in the derivation of Eqn. (0). Fix also a field S" F on 



S, and a positive number e. We wish to show existence of a smooth in 
l^, with = that is an "e-approximate" solution of (p3|), i.e., that is 
such that the square root of the integral over V of the square of the left side 
is less than or equal to e. We may set 0o = (by replacing by + in 
(|93D, where 0|5 = 0o, and then absorbing into ja the extra term, pS/ m4>^ i 
thus created). Then, to show existence of an e-approximate solution of (|93|) , 
it suffices to show that fields of the form k^^ m'P^ ) with smooth and 
vanishing on S", are L^-dense in V . But for this, in turn, it suffices to show: 
Given any square-integrable field '0° on V , such that 



V 



rkorpVn.<p^ = (94) 



for every smooth vanishing on 5, we must have if)"" = 0. So, fix such a 
-0°. Were this ijj smooth, then we could proceed as follows. Integrate the 
left side of Eqn. (0) by parts. First choosing to have support in the 
interior of V, we obtain (from the volume integral) that V m{ka^ = Oj 
and then choosing to be nonzero on St„, we obtain (from the surface term) 
that ip^\sto ~ 0- these two together imply that -0" = 0. Indeed, setting 
^L'"/3 = koTp , L = i^mka'^p)^", j'a = 0, = V^, 0' = in (0), we obtain 



l^a 'V ^mabac 

' Sto 

+4(Tr||(Vl)^/'e-*/^|| +4(Tr||V(^^)e-*/^||, (95) 

while, for r sufficiently small, the second term on the right is less than or 
equal to the left side. But, unfortunately, -0° need not be smooth, but only 
square-integrable. We therefore proceed as follows. Let hhe a, nonnegative, 
smooth, symmetric, two-point function on M, such that h{x,y) vanishes for 
X and y sufficiently separated, and, for each fixed y, the integral of h{x, y) 
over M has value one. Set '0(x) = /^^ h{x,y)ip{y), where the integration 
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variable is ?/ G M. This is a smooth approximation to ij). Now take a 
sequence of such /I's such that "/i, together with its derivative, approaches a 
delta-function and its derivative" , in the sense that 



hix,y)f{y) ^ fix), V h{x,y)f{y) ^Vfix), (96) 

My J My 

for every smooth function / on M. Then the left side of (p^), with re- 
placed by tp, approaches zero, and so, it can be shown, each of the first and 
third terms on the right in (|95|) , with ip replaced by ip, approaches zero. It 
now follows from ( pS]) (again, choosing r sufficiently small) that | |?/'e~*/'^| | 
approaches zero, and so that = 0. Thus, we have shown existence of an 
e-approximate solution of Eqn. (p3D, in a certain region V, with given initial 
data. 

We now adopt an iterative procedure. Begin with any smooth field 
01 on V, satisfying the initial condition, ^il^ = 0o- Choose e2 > 0, set 
ka"^/3 = ka^fi{x, 01 (x)), ja = ja{x, 4>i{x)) in (pSj), and find an e2-approximate 
solution, 02, of that equation satisfying the initial condition. Then choose 
£3 > 0, set k = /c(02), j = j(02) in (P5|), and find an ea-approximate solu- 
tion, 03, of that equation satisfying the initial condition. Continuing in this 
way, with — 0, we obtain a sequence of fields, 0i, 02, ■ ■ ■. Each of these 
satisfies the initial condition, and each approximately satisfies (0) (better, 
as i increases) with its predecessor as background. We wish to show that 
the 0j converge to the desired solution. To this end, set = 0j, 0' = 0j_i, 
k = A;(0i_i), k = A;(0i_2), J = j(0i-i), / = j(0i-2) in (0), to obtain 

II0.-0.-1II < are*°/^[2||(A;(0,_i)-A;(0,_2))V(0, + 0,_i)|| 
+4||j(0.-i)-j(0.-2)||+4(e, + e,_i)] 
< are*°/^[a||0,_i-0,_2||+4(e, + e,_i)], (97) 

where we used e~^°l^ < e"*/"^ < 1 in the first step; and set a = 2 lub||^|lub 
|V(0i + 0t-i)| + 4 lub|||| in the second. Suppose for a moment that 0j and 
V0i were uniformly bounded (i.e., there is a single constant that bounds all 
the 0i throughout V, and similarly for V0j). Then a (via (|9^) ) and a (above) 
would remain bounded as i increases, and so, by choosing r, to, and the 
sufficiently small[^ in (0), we would guarantee convergence of J2 I |0i — 0i-i| I, 
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It is here that, by having to choose to sufRciently small, we restrict the size of the 
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and so convergence of the 0j to some field (f) on V. But, unfortunately, we 
cannot, at this stage, guarantee uniform boundedness of even the 0i — much 
less of their derivatives. Indeed, all we control is a certain average of the (pi, 
represented by the left side of (|97D. We therefore proceed as follows. Taking 
the x-derivative of Eqn. (|89D , and introducing a new field to represent 

"V0" , we obtain 

ka'^pVmi(p/) + (pJiVakoT'p + 0a^^fc„™^) + V aja + 0a^^Ja = 0, (98) 

The combination of (|89D and (|98D is a quasilinear, first-order system^ of par- 



tial differential equations for the fields 0", (po!^. Furthermore, it inherits from 
(|89D its hyperbolicity (since the coefficient of the derivative-term in (|98D is 
the same k as in (^)). Now continue in this way, taking successive deriva- 
tives of (^9]), introducing successive fields, (pa---c°', to represent the higher 
derivatives of 0°, and obtaining successively larger hyperbolic systems. Con- 
sider the system that results after taking the fourth derivative — so the fields 
are now 0°,0a°, 0a6", 0a6c", (pabcd^ and the equations (p9| ) and its first four 
derivatives — and apply to this entire system the iterative procedure above. 
Then the left side of (|97| ) will include a term ||0j abed" — (pi-i ahccf'W- We now 
apply a Sobolev inequality^, which asserts that ||VVVV0|| provides a uni- 
form bound on and its first two0 derivatives. From this, combined with 
(1^, it can be shown that 0i", 0ia°, (piah^, 4>iabc", (piahcd^ all converge, to 
some fields 0", 0a", 0afe", 0a6c", (pabcd"- The 0" that results from this proce- 
dure is the desired solution of (pQ]). To show that this 0" is smooth, take still 
higher derivatives of (^), and proceed as above, introducing as new fields 
still higher derivatives of 0, and applying the interative procedure above to 
the resulting hyperbolic system. One must check that, in the demonstra- 
tion of existence for this succession of hyperbolic systems, the number to 

neighborhood U for part (a) of the theorem. (In the large-t region of V , the could fail 
to converge, indicating that there the final solution would become singular.) This is the 
starting point for deriving a strengthening of part a guaranteeing a minimum size to its 
neighborhood U. 

^^For initial data on 5, we take, for </>", the given 0O: and, for 0^", "Va(/>"" , as computed 
from 00 and Eqn. (|9|) evaluated on S. 

^"^R. A. Adams, Sobolev Spaces (Academic Press, NY, 1975). 

^^We need, at this point, a uniform bound on the second derivative of <j) because there 
appears in our equations, by this point, a term "(VV0)(VVV0)", and we have on VVV0 
only an bound. 
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(which governs the size of the region V) can remain bounded away from 
zero. There results convergence of cpia---" = Va---Vc0i", and so smooth- 
ness of 0" = hm 
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